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ABSTRACT 


\ 

Most  aerospace  structural  materials  exhibit  some  degree  of  aniso¬ 
tropic  strain  hardening.  During  the  past  few  years,  several  methods 
have  appeared  in  the  literature  for  introducing  inelastic  isotropic 
material  behavior  effects  into  existing  matrix  analysis  procedures  using 
the  incremental  theory  of  plasticity.  A  review  is  presented  of  these 
methods  and  a  step-by-step  routine  known  as  the  "constant  strain"  method 
is  selected  for  the  development  of  an  anisotropic  inelastic  procedure. 

«* 

A  simple  truss  with  one  redundant  is  used  to  indicate  the  basic 
ideas  of  the  approach.  Then  the  procedure  is  generalized  to  the  more 
important  case  of  biaxially  stressed  structures.  Nodal  stresses  are 
evaluated  step-wise  for  increasing  load  through  the  use  of  an  influence 
coefficient  equation.  '  The  inelastic  (plastic  and  creep)  strains  at  one 
load  level  are  used  as  initial  strains  at  the  subsequent  level  to  account 
for  nonlinear  effects.  The  anisotropic  behavior  is  considered  by  using 
a  proposed  extension  of  Hu's  strain  hardening  theory. 

Several  analyses  of  an  aluminum  alloy  (202h-Th)  shear  lag  structure, 
which  has  been  tested  previously  for  the  Air  Force,  are  carried  out, 
first  assuming  isotropic  and  then  anisotropic  material  properties.  The 
correlation  between  test  results  and  those  predicted  by  isotropic  theory 
is  reasonably  good.  The  anisotropic  analysis  gives  predicted  results 
which  are  in  slightly  more  consistent  agreement  with  the  test  data. 

The  procedure  is  also  modified  to  give  an  isotropic  deformation 
theory  solution,  which  produces  numerical  results  in  a  much  shorter  com-, 
puter  time  than  required  for  the  incremental  theory  solution.  In  the 
case  of  the  shear  lag  structure  investigated,  the  results  by  the  two 
theories  are  in  very  close  agreement. 

Creep  test  results  of  an  1100-F  aluminum  shear  lag  structure  are 
also  available.  An  analysis  of  this  structure  by  the  proposed  incre¬ 
mental  method  is  carried  out  and  its  predictions  too  are  in  reasonably 
good  agreement  with  the  test  data.  The  1100-F  material  is  very  nearly 
isotropic  and  no  testing  of  structures  exhibiting  anisotropic  creep  is 
known  to  have  been  performed.  Hence  the  anisotropic  creep  capability 
of  the  proposed  method  cannot  be  checked  out  against  tests  at  this  time. 

A  sample  calculation  is  nevertheless  carried  out  for  a  hypothetical 
material  having  this  characteristic. 

The  approach  presented,  which  is  simple  in  concept  and  execution, 
is  found  to  be  a  reasonably  good  phenomenological  model  of  an  exceedingly 
complex  physical  problem.  The  accompanying  digital  computer  program  is 
believed  to  be  very  versatile,  and  well  suited  for  the  inclusion  of  any 
other  types  of  material  nonlinearity  that  may  be  of  interest. 
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SECTION  I 
INTRODUCTION 


With  the  introduction  of  high-speed,  large-capacity,  digital  computers, 
a  number  of  investigators  (References  1,4,7)  have  adapted  essentially 
linear  matrix  analysis  methods  to  the  solution  vf  redundant  structures 
where  nonlinear  material  plasticity  and  creep  properties  are  considered. 
These  methods  have  been  either  iterative  or  non-iterative  step-by-step 
numerical  procedures. 

The  present  work  is  an  effort  to  explore,  revise,  and  extend  the 
matrix  method  of  analysis  in  order  to  apply  it  to  a  range  of  practical 
aerospace  structural  problems  exhibiting  inelastic  isotropic  or  inelastic 
anisotropic  material  behavior.  Thus,  the  intent  is  to  concentrate  upon 
methods  which  are  able  to  predict  inelastic  strain  distributions  in 
irregular  idealized  structures  in  a  biaxial  stress  state  where  materials 
exhibit  strain  hardening  and  creep  properties  representative  of  those  em¬ 
ployed  in  aerospace  construction. 

A  discussion  of  complete  load  reversal,  although  desirable  for 
plastic  fatigue  studies,  is  not  included  because  the  theoretical  founda¬ 
tion  for  such  a  procedure  is  apparently  not  yet  fully  developed.  The 
present  formulation  of  elastic,  plastic  and  creep  loading,  followed  by 
elastic  unloading,  while  restrictive,  is  nevertheless  of  practical  inter¬ 
est. 

Of  the  proposed  analytical  methods,  one  by  Denke  (Reference  l)  has 
developed  naturally  frcm  the  matrix  force  method  of  analysis  and  consists 
of  including  nonlinear  plastic  and  creep  terms  in  the  equations  for  the 
gaps  at  those  cuts  which  are  required  to  make  the  structure  statically 
determinate.  The  redundant  forces,  required  to  close  the  gaps  and  make 
the  structure  continuous,  are  obtained  by  solving  these  nonlinear  equations 
by  a  Newton-Raphson  procedure. 

A  second  method  reported  by  Kobayashi  and  Weikel  (Reference  2)  has 
been  developed  from  the  direct  stiffness  method.  Here,  forces  occurring 
as  a  consequence  of  the  inelastic  effects  are  included  in  the  nodal  force 
equations.  The  nodal  displacements  (or  displacement  rates,  where  creep 
is  considered)  obtained  by  solving  these  equations  impose  equilibrium 
at  the  nodes  under  the  action  of  internal  loads,  surface  tractions,  and 
prescribed  displacements . 

When  the  inelastic  effects  are  accounted  for  by  a  flow  (or  incremental) 
theory,  the  deformation  is  an  accumulation  of  increments  each  governed  by 

the  prevailing  stress.  Thus,  when  either  of  the 'methods  just  described 
is  used,  one  must  obtain  (for  the  governing  equations)  a  series  of  solu¬ 
tions,  with  one  solution  corresponding  to  each  load  increment.  This  re¬ 
quires  a  considerable  amount  of  computer  time  even  for  medium-sized 
analyses. 
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The  approach  recommended  here,  in  addition  to  being  conceptually 
simple,  does  not  require  repeated  matrix  inversions.  It  was  developed 
fran  that  proposed  oy  R.  Gallagher,  J.  Padlog  and  others  at  Bell  Aero- 
'system  (References  Input  data,  as  in  the  case  with  these  other 

nonlinear  analyses,  are  generated  by  an  elastic  analysis;  however,  for 
this  approach,  either  the  matrix  force  method  or  direct  stiffness  method 
can  be  used.  The  problem,  formulated  in  terms  of  standard  influence 
coefficients  for  applied  load  and  initial  strain,  is  reduced  from  a  non¬ 
linear  to  a  linear  one  by  using  those  strains  obtained  at  the  previous 
load  level  to  approximate  the  current  inelastic  strains. 

Development  of  the  anisotropic  analysis  is  based  on  an  extension 
of  the  proposed  anisotropic  theory  of  Hu,  Reference  11.  The  constant- 
anisotropic-coefficient  assumption  of  Hu  is  replaced  by  one  in  which 
the  coefficients  are  allowed  to  vary  with  the  level  of  stress.  The 
formulation  is  then  a  simple  modification  of  the  isotropic  procedure. 

It  is  also  shown  that  anisotropic  creep  can  be  included  in  a  manner 
similar  to  the  isotropic  creep. 


SECTION  II 


INELASTIC  MATRIX  METHODS 


A.  Formulation 

An  inelastic  structural  analysis  can  be  carried  out  in  two  steps. 

The  first  is  the  standard  elastic  solution  where  internal  stresses  and 
accompanying  strains  are  related  through  Hooke's  law.  The  second  step, 
the  modification  of  this  elastic  system  to  include  inelastic  strains, 
is  analogous  to  a  procedure  for  including  superimposed  thermal  strains. 
The  inelastic  strains  are  defined  as  the  differences  between  the  total 
strains  and  the  elastic  strains  and  are  generally  functions  of  the  final 
stresses,  not  those  of  the  linear,  elastic  state. 

The  simple,  pin-jointed,  single-redundant  truss,  pictured  in  Figure 
(l),  illustrates  the  basic  notions  more  clearly.  All  bars  are  consid¬ 
ered  to  be  elastic,  except  the  vertical  diagonal  which  can  became  plastic 
and  has  a  stress-strain  relation  represented  by  the  curve  shown  in  Fig¬ 
ure  (2).  The  applied  load  P  is  large  enough  to  cause  member  3  to  become 
plastic . 

A  solution  might  be  obtained  by  first  simply  ignoring  plastic  strain 
in  member  3  and  assuming  all  members  elastic.  The  resulting  stress 
would  then  be  of  magnitude  a.  The  actual  stress  for  member  3,  of  lower 
magnitude  due  to  plastic  yielding,  is  designated  and  the  associated 

inelastic  strain  in  Figure  (2).  These  stresses  and  strains  are 

related  in  the  following  equation. 


a  + 


(1) 


where: 

(k) 

a  is  the  actual  final  stress 

a  is  the  elastic  stress 

(k) 

e  the  inelastic  strain 

T  the  redundant  elastic  stress  for  a  unit  value  of 
initial  strain 


To  be  more  specific  T  is  equal  to  the  redundant  stress  in  member  3 
corresponding  to  a  unit  initial  strain  in  member  3.  Note  that  T  must 
be  negative  to  cause  a  reduction  of  stress  in  the  diagonal  member. 

An  important  feature  to  be  observed  about  Equation  1  is  that, 
since  the  inelastic  strain  e(k)  is  a  nonlinear  function  of  the  final 
stress  oM,  this  equation  is  really  a  nonlinear  relation  to  be  solved 
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(k) 

for  a  ' .  This  characteristic  will  always  he  present  in  the  analyses 
to  he  discussed  in  this  report. 

Equation  1  can  be  generalized  to  provide  stresses  in  all  members 
of  the  structure  and  provide  for  inelastic  strains  in  all  members  for 
a  variety  of  applied  loads.  This  basic  influence  coefficient  equation* 
is  as  follows: 


*  truv]f^) 


(2) 


The  ct  'a  are  the  ordinary  stress  components  at  the  various  node 

points  of  the  structure.  An  element  of  [T  ]  gives  the  u  stress  in 

bm  th 

the  linear  redundant  structure  due  to  a  unit  in  applied  load,  and 
{Pm }  represents  the  actual  applied  loads.  Also,  an  element  of  [I*  ] 

gives  the  uth  stress  component  in  the  linear  redundant  structure  due 

to  a  unit  initial  strain  at  the  v  stress  location  in  the  unloaded, 
statically  determinate  structure.  Finally,  an  element  of  [e  }  repre- 

sents  the  actual  initial  strain  at  the  v6  stress  location. 

The  problem  is  now  reduced  in  essence  to  the  determination  of  the 
inelastic  strains  to  use  as  the  initial  strains  {ev}  in  Equation  2. 

For  a  structure  in  a  load-temperature-time  environment,  this  task  can 
be  rather  formidable,  because  { is  a  function  of  local  temperature 

and  time,  as  well  as  the  local  stress  history. 

B.  Example  Problem 

Continuing  our  study  of  the  simple  truss  example  of  Figure  (l),  we 
now  allow  all  members  to  go  plastic.  The  exact  results  for  the  defor¬ 
mation  and  stresses  in  the  truss,  for  nonlinear  properties,  are  easily 
obtained  by  direct  numerical  solution  of  the  equations,  and  hence  will 
be  used  without  development. 

The  step-by-step  finite  element  method  for  determining  the  stresses 
{au}  and  strains  {e^.}  involves  the  use  of  Equation  2  and  a  nonlinear 

stress-strain  relation.  This  relation  will  be  assumed  to  be  a  piece- 
wise  linear  approximation  of  the  Ramberg-Osgood  stress-strain  relation. 


•The  derivation  of  matrices  [r  ]  and  [r  ]  is  given  in  Appendix  VI. 
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vhere  e  denotes  the  Inelastic  (or  plastic)  strains  and  is  given  by 


and  where 

£  =  total  strain 
a  =  member  stress 
E  =  Young's  modulus 

a  =  reference  stress  (stress  at  secant  modulus  of  0.7E) 
o 

0  =  nonlinear  parameter 

The  first  step,  in  applying  the  finite  element  method,  is  to  obtain 

the  influence  coefficient  matrices  [T  ]  and  [r  ]  for  the  linear,  re- 

um  uv 

dundant  structure.  This  requires  specification  of  the  geometry  of  the 
structure  and  the  (linear)  material  properties  of  the  individual  struc¬ 
tural  elements.  The  geanetry  of  the  example  truss  problem  is  given  in 
Figure  (l).  Th  j  material  is  assumed  to  be  an  aluminum  alloy  with  the 
following  constants. 


E  =  10^  psi,  aQ  =  105  psi,  0  =  10 


Hie  approximate  stress -strain  curve  used  matches  the  Rsmberg-Osgood  values 
at  2000  psi  stress  intervals. 

In  this  case,  we  find 


CW  - 


rr  ]  =  io' 

L  uvJ 


[■ 


klk 

.207 

.207 

586 

-  .293 

-  .293 

586 

-  .293 

-  .293 

for  the  case  of  the  single  applied  load  P 


C.  Step-by-Step  Methods 

References  3  and  1*  present  what  appears  to  be  the  simplest  possible 
approach  to  this  problem  fran  a  computational  standpoint.*  A  non  iter¬ 
ative  step-by-step  calculation  is  performed  in  which  all  quantities  in¬ 
cluding  the  initial  strains  ^  are  incremented  and  then  assumed  to  remain 

constant  in  the  ensuing  load  interval.  The  inherent  difficulty  in  this 
approach  is  to  establish  the  connection  between  successive  steps.  Two 
methods,  both  of  which  involve  the  initial  strains  from  the  prior  step 
to  predict  quantities  in  the  current  step  are  suggested.  It  is  antici¬ 
pated  that  by  controlling  the  size  of  the  interval  one  may  achieve  any 
degree  of  accuracy. 

The  development  herein  is  discussed  only  in  the  detail  necessary  to 
analyze  the  redundant  truss  of  Figure  (l)  being  loaded  for  the  first  time 
into  the  plastic  range.  Generalization  to  biaxial  plasticity  and  creep  phe 
nanena  are  discussed  in  the  succeeding  sections. 

The  step-by-step  procedure  for  solving  the  problem  is  introduced  by 
rewriting  Equation  2  in  the  form 

k00}  ■  ♦  Eruv4v<k'1)}  <3) 

where 

k  is  the  cycle  designation 


This  can  be  regarded  as  the  fundamental  equation  for  the  non-iterative, 
step-by-step  methods.  The  idea  in  formulating  this  equation,  as  indicated 
by  the  cycle  designating  superscript  is  that  the  initial  strains  of  the 
previous  cycle  can  be  used  to  approximate  the  initial  strains  of  the 
current  cycle.  The  strains  of  the  previous  cycle  may  be  incorporated  in 
several  ways,  two  of  which  constitute  the  constant  stress  and  the  constant 
strain  methods  of  analysis. 

D.  Constant  Stress  Method 

As  indicated  previously,  in  the  step-by-step  procedure  considered 
here,  one  enters  the  k^*1  cycle  with  applied  loads  jp^^)  |  and  initial 
strains  latter  evaluated  during  the  preceding  cycle. 


•These  methods  make  use  of  devices  previously  used  by  others  to  solve 
inelastic  problems;  for  example,  S.  S.  Hanson  at  the  Lewis  Research 
Laboratory,  NASA,  Cleveland,  Ohio,  has  previously  carried  out  inelastic 
analyses  of  turbine  discs  involving  somewhat  similar  techniques 
(Reference  5)* 
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The  first  operation  of  the  current  cycle  is  to  determine  frco 

Equation  3  by  direct  substitution.  The  second  operation  is  a  determination 
of 

{6v(k)}  for  use  in  the  next  cycle.  The  constant  stress  method  does 

this  in  the  most  obvious  way,  by  reading  from  the  given  stress-strain 

curve  the  plastic  strains  |ev^  j  corresponding  to  the  cr^^'s  (the 

reason  for  the  name  "constant  stress"  is  thus  apparent).  The  operation 
is  indicated  schematically  in  Figure  (2). 


Ml 


The  results  of  the  application  of  this  method  to  the  example  truss 
problem  are  shown  in  Figure  (4),  where  the  stress  in  the  vertical  member 
(Bar  #3)  has  been  plotted  versus  load.  These  results  display  a  striking 
defect  of  the  method  due  to  the  development  of  a  sudden  and  catastrophic 
divergence,  whose  onset  depends  upon  step  size.  This  dependence  is  such 
that  any  attempt  to  improve  accura  jy  reducing  step  size  only  hastens 
the  occurrence  of  divergence.  An  c  ..sanation  of  this  behavior  is  given 
in  Appendix  III.  Because  of  this  defect,  the  constant  stress  method  in 
this  form  must  be  eliminated  frcm  consideration  as  an  acceptable  method 
for  general  use. 

E.  Constant  Strain  Method 

The  first  operation  of  the  constant  strain  method  is  exactly  the 
same  as  the  first  operation  of  the  constant  stress  method;  a  W  is 

evaluated  by  direct  substitution  in  Equation  3*  Thereafter,  one  de¬ 
termines  «  ^  for  use  in  the  next  cycle  as  follows.  Referring  to 
Figure  (3),  for  each  member,  point  A  is  determined  with  stress-strain 
coordinates  <j  and  av  M/B  +  A  relaxed  stress  a  *00  is 

now  calculated  with  the  same  total  strain,  corresponding  to  point  B  on 
the  given  stress-strain  curve.  Note  that  here  the  total  strain,  rather 
than  stress,  remains  unchanged— hence  the  name  "constant  strain"  method. 

The  required  initial  strain  is  the  inelastic  strain  cM  correspond¬ 

ing  to  the  relaxed  stress,  as  indicated  on  Figure  (3). 

The  results  of  applying  the  constant  strain  method  to  the  truss 
problem,  for  the  three  step  sizes  5000,  500  and  50  lb.,  are  shown  in 
Figure  (5).  The  accuracy,  for  a  given  step  size,  is  not  as  good  as 
that  of  the  constant  stress  method,  but  the  analysis  is  now  free  of  any 
instability.  The  constant  strain  method  is  therefore  selected  for 
further  use  herein.  The  discussion  of  the  step  size  and  of  a  method  of 
monitoring  it  is  left  for  a  later  section. 


SECTION  in 


ISOTROPIC  ELASTIC-PLASTIC  ANALYSIS 


A.  Biaxial  Theory 

Having  presented  the  simple  truss  example  of  the  step-by-step  pro¬ 
cedure,  we  proceed  now  to  the  case  of  more  practical  interest  —  a  bl- 
axially  stressed  structure.  The  new  procedure  is  identical  to  the  one 
already  discussed  for  the  simple  truss  with  one  exception.  Because  of 
the  biaxial  stress  we  can  no  longer  work  directly  frcm  the  stress-strain 
curve  to  obtain  the  plastic  strains  for  use  in  Equation  3J  instead,  we 
must  employ  the  well-known  concept  of  an  "effective"  stress-strain  re¬ 
lationship  in  conjunction  with  a  von  Mises  type  yield  condition  and  the 
associated  incremental  flow  relations. 

The  biaxial  theory  is  described  by  a  summary  of  the  steps  to  be 
used  as  a  guide  for  a  detailed  description  which  follows.  The  constant 
strain  method  used  here  is  a  step-by-step  procedure  which,  after  incre¬ 
menting  the  applied  load,  can  be  applied  in  four  parts: 

1.  Obtain  the  stress  components  at  each  node  using  the  basic 
Equation  3  by  assuming  the  initial  strains  fron,  tne  previous 
load  level. 

2.  Using  these  stresses,  calculate  an  effective  stress  at  each 
node. 

3«  Assume  that  the  effective  stress-strain  relation  for  the 

material,  modified  by  including  the  elastic  strain,  corresponds 
to  data  measured  in  a  simple  uniaxial  tension  test.  Using  this, 
calculate  the  effective  strain  corresponding  to  the  effective 
stress. 

4.  Using  the  incremental  flow  relations,  determine  the  inelastic 
strain  increments.  The  proportionality  constant  in  these 

equations  is  the  ratio  of  the  effective  strain  increment  to  K 

the  effective  stress. 

At  this  point  in  the  calculation,  the  applied  load  can  be  incremented 
again  and  the  cycle  repeated. 

When  calculating  the  ordinary  stresses  for  the  kth  load 

level  using  Equation  3  (Step  l),  it  is  convenient  to  re- identify  these 
stresses  by  means  of  a  new  subscript  N,  as  follows: 


8 


r  •  -\ 


3H-2 

,00 

3K-1 

,00 

r3H 


(4) 


The  stresses  are  thus  arranged  in  groups  of  three  components  (two  normal 
and  one  shear)  at  each  node  point  N. 

(k) 

We  now  calculate  the  corresponding  effective  stresses  '  for  each 
of  the  nodes  from  the  von  Mises  type  formula  (Step  2) 


Note  that  by  this  definition 


must  be  positive  and  is  proportional 


to  the  octahedral  shear  stress.  This  formula  together  with  the  stress- 
strain  data  constitutes  the  strain  hardening  criterion. 


We  now  go  to  the  tensile  stress-strain  curve  (which  is  also  the  9  vs 
o/E  +  e(p)  curve)  for  the  material  of  interest  and,  using  the  constant 
strain  method,  read  from  it  the  corresponding  effective  plastic  strain 
;(*)'*>.  This  operation  (Step  3)  ia  identical  to  that  previously  de¬ 
scribed  for  the  uniaxial  case  on  page  7. 


In  accordance  with  the  flow  theory  of  plasticity,  the  increment  in 

the  effective  strain  Ai(p)<k)  over  that  of  the  preceding  interval  must 

be  calculated  (Step  4).  The  increment  will  be  either  positive  or  zero, 
depending  upon  whether  plastic  loading  or  elastic  unloading  (or  reloading) 
is  taking  place.  Thus, 


4I(p41)  -  -«(p)|jlt)  -  J(p)£k-1)  (6a) 


when 


is  greater  than  any  previous  5^  (inelastic  strain  increasing) 
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1  *  * 


(6b) 


Ae(p4k)  =  0 

when  5^  is  smaller  than  a  previous  &N  (elastic  unloadir  jv  reloading 
inelastic  strain  constant) 

The  increments  in  the  ordinary  plastic  strain  components  may  now 
be  obtained  using  a  Prandtl-Reuss  incremental  relationship. 


The  total,  ordinary  plastic  strain  components  are  obtained  by 
addition. 


These  components  together  with  the  new  applied  loads  P^kf1^ 

may  be  substituted  in  Equation  3  bo  obtain  in  the  next  load 

cycle. 


10 


B.  Determination  of  Calculation  Step  Size 


It  should  be  noted  that,  according  to  the  constant  strain  method, 
every  predicted  value  of  effective  strain  .  ?<4t!  together  with  its 
accompanying  value  of  effective  stress  aj^)  constitutes  an  approxi¬ 
mation  to  a  point  on  the  actual  effective  stress-strain  curve.  The 
excellence  of  the  approximation  is  directly  related  to  the  loading  in¬ 
crement,  as  is  shown  in  the  truss  results  Figure  ( 5 ) *  Thus  it  is  only 
necessary  to  monitor  this  agreement  for  one  or  more  of  the  critically 
loaded  nodes  to  determine  whether  the  step  size  is  satisfactory.  This 
is  illustrated  below  in  connection  with  the  shear  lag  structure  investi¬ 
gation. 

C.  Description  of  Shear  Lag  Structure 

Several  very  useful  tests  have  been  performed  for  the  Air  Force 
upon  shear  lag  structures  (Reference  7) •  The  structure,  loaded  as  shown 
in  Figure  (6),  is  an  integrally  machined  part  of  2024-T4  aluminum  alloy 
stiffened  along  the  loading  (y)  axis.  The  stiffener  is  tapered  in 
thickness  from  each  end  towards  the  center  of  the  structure. 

This  structure  was  chosen  originally  because  it  is  simple  to  work 
with  and  .:ell  adapted  to  analysis  by  both  matrix  methods  when  appropri¬ 
ate  idealizations  are  employed.  When  tension  forces  are  applied  to  the 
ends  of  the  stiffener,  high  stress  gradients  are  induced  in  a  manner 
analogous  to  those  encountered  in  aircraft  structures. 

The  material  properties  essential  to  this  analysis  were  obtained 
from  tension  tests  reported  in  Reference  7»  These  tests  were  performed 
on  coupons,  machined  from  the  parent  plate,  in  the  longitudinal  or 
x-direction  and  the  transverse  or  y-direction  of  Figure  (6).  The  data 
resulting  from  these  tests.  Figure  (ll),  indicate  the  presence  of  a 
considerable  degree  of  anisotropy.  In  the  present  study,  three  piece- 
wise  linear  representations  of  stress-strain  curves  were  fitted  to  these 
points;  two,  R01  and  R02,  in  Figure  (ll),  are  equivalent  to  Ramberg-Osgood 
curves  used  in  Reference  7;  the  third,  R02M,  is  a  Grumman  modification. 

The  modulus  of  elasticity  of  all  the  curves  is  taken  as  10.3  x  1  (P  psi. 
Note  also  that  the  maximum  strains  recorded  are  of  the  order  of  0.010 
in/in,  whereas  the  maximum  strains  reached  in  the  shear  lag  tests  are 
around  0.020  in/in.  Thus  there  is  some  doubt  as  to  whether  the  idealized 
curve  correctly  represents  the  test  material  in  this  high  strain  region. 


The  locations  of  the  strain  gages  for  the  test  of  the  stiffened 
plate  are  shown  in  Figure  (13) •  The  plate  was  loaded  by  applying  tension 
to  the  stiffener  in  steps  of  1000  pounds  to  6000  pounds,  gage  readings 
being  taken  at  each  step.  It  was  then  unloaded  in  steps  of  1500  pounds 
to  zero,  and  finally  progressively  loaded  to  failure.  Buckling  occurred 
at  a  load  of  23,000  pounds  and  fracture  at  25.800  pounds.  Data  from  this 
test  are  plotted  on  Figures  (l4)  through  (22). 
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D.  Elastic  Shear  Lag  Structure  Analysis 


The  idealizations  of  the  upper  right  guandrant  of  the  shear  lag 
structure  for  a  direct  stiffness  and  a  force  method  analysis  are  shown 
in  Figures  (7)  and  (8)  respectively.  A  typical  element  of  the  force 
method  analysis,  Figure  (9),  consists  of  conventional  bars  and  rectangu¬ 
lar  shear  panels.  Many  previous  idealizations  have  emitted  the  Poisson's 
Ratio  effect.  The  present  idealization,  however,  incorporates  this 
effect  in  the  manner  described  in  Appendix  II. 

The  basic  element  of  the  stiffness  approach  consists  of  a  cluster 
of  four  "Turner  triangles"  (References  8,9)  to  form  a  rectangle  as  shown  in 
Figure  (10).  The  manner  of  obtaining  the  stresses  is  discussed  in 
Appendix  I. 

An  elastic  analysis  under  a  unit  applied  load  was  performed  by  both 
the  force  and  the  stiffness  methods.  These  results  are  compared  in 
Appendix  I. 

The  inelastic  analysis  can  be  made  using  either  of  the  two  approaches 
(stiffness  or  force  method).  For  the  present  investigation,  only  the  force 
method  is  used. 

E.  Flow  Theory  Shear  Lag  Structure  Analysis 

The  inelastic  analysis  was  carried  out  using  the  Fortran  2  program 
listed  in  Appendix  V.  This  program  is  capable  of  carrying  out  isotropic 
or  anisotropic,  plastic  or  creep  flow  theory  analyses.  The  flow  charts 
and  instructions  for  preparation  and  submission  of  data  are  also  included 
in  Appendix  V. 

Before  comparing  the  analytical  and  test  results,  let  vis  look  at  a 
plot  of  the  tensile  stress-strain  curve,  R034,  and  compare  It  with  the 
predicted  effective  stress-strain  relationship  for  various  step  sizes. 

Such  a  comparison  is  found  in  Figure  (12)  for  the  node  corresponding  to 
the  center  of  the  specimen  vhich  is  the  point  of  highest  strain  in  the 
structure .  It  can  be  seen  there  that  for  a  step  size  of  AP  =  500  lh., 
the  agreement  Is  rather  poor.  For  AP  =  50  lb.  the  agreement  is  much 
better,  while  for  AP  =  5  lb . ,  the  predicted  value  lies  directly  on  the 
stress-strain  curve.  The  TEH  7094  computer  time  for  this  best  result  and 
a  maximum  load  of  P  =  16,760  lb.  is  approximately  20  minutes. 

The  predicted  strain  distributions  are  shown  in  Figures  (l4)  through 
(2l),  together  with  the  corresponding  test  values,  along  the  two  strain 
gage  lines.  In  these  plotB,  the  calculated  results  are  linearly  inter¬ 
polated  values  between  node  points.  Figures  (l4)  and  (15)  give  elastic 
results;  the  agreement  with  test  data  is  seen  to  be  rather  good,  giving 
the  necessary  confidence  in  the  accuracy  of  the  basic  Influence  co¬ 
efficient  matrices.  It  is  observed  also  that  the  specimen  achieves  its 
basic  purpose  of  displaying  a  pronounced  shear  lag  effect  with  the 
highest  strain  occurring  at  the  central  node,  as  expected. 
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As  the  applied  load,  increases  through  7070  lbs.,  the  strains  at  the 
central  node  become  plastic.  During  the  tests,  the  strain  gages 
continued  to  function  through  an  applied  load  level  of  14,600  lb.;  beyond 
this  point  the  x-gage  failed  to  record.  The  y-gage  failed  also  above  a 
load  level  of  16,760  lb.,  which,  consequently,  is  the  highest  level  con¬ 
sidered  in  the  comparison  of  test  and  analysis  even  though  the  plate  did 
not  buckle  until  P  =  23,000  lb. 

The  proportional  limit  for  the  tea'  specimen  material  occurs  at  a 
strain  of  approximately  .004  in/in,  as  shown  in  Figure  (ll).  Bearing  this 
in  mind  during  an  examination  o^  Figures  (l6)  through  (2l),  it  is  seen  that 
plastic  behavior  is  primarily  confined  to  a  fairly  small  region  around 
the  central  node,  and  we  thus  have  a  case  of  contained  plasticity.  The 
analysis  predicts  a  very  pronounced  strain  redistribution  extending  some¬ 
what  beyond  this  region.  This  is  indicated  by  a  comparison  of  the  plastic 
results  with  extrapolated  elastic  results  shown  as  dotted  curves  on  Fig¬ 
ures  (20)  and  (2l). 

As  for  agreement  with  test  values,  the  analysis  substantially  undfer- 
predicts  the  strain  gage  readings  where  plasticity  is  most  pronounced. 

Since  the  elastic  results  agree  so  much  better,  one  must  assume  that  the 
difficulty  lies  somewhere  in  the  plasticity  part  of  the  correlation.  It 
was  mentioned  previously  that  the  idealization  of  the  stress-strain  curve 
R02M  of  Figure  (ll)  was  open  to  question  at  the  high  strain  end  because 
of  the  absence  of  test  points.  Accordingly,  an  additional  run  was  made 
for  a  revised  curve  extending  horizontally  beyond  the  last  indicated  test 
point.  The  plastic  strains  at  the  critical  central  node  are  increased 
approximately  10$  by  so  doing.  This  represents  an  appreciable  closing 
of  the  gap,  but  the  gap  nevertheless  remains.  N 

Calciilations  were  also  made  based  upon  the  previously  mentioned  stress- 
strain  curves  R01  and  R02  of  Figure  (ll).  The  results  for  the  central  node 
are  shown  on  Figure  (22).  As  might  be  expected,  they  depart  appreciably 
fraa  the  R02M  predictions. 

F.  Deformation  Theory  Analysis 

Solutions  by  a  deformation  theory  have  traditionally  been  considered 
to  be  more  easily  obtainable  than  flow  theory  solutions.  This  is,  of 
course,  because  only  the  stresses  at  the  final  applied  load  level  need 
be  considered,  rather  than  the  stress  histories  developed  during  loading. 

It  is  therefore  of  interest  to  determine  whether  similar  benefits  are 
attainable  in  the  case  of  the  finite  element  analyses  currently  being 
considered. 

Once  again,  a  solution  of  Equation  3  is  required,  this  time  such  that 
the  initial  strains  e  satisfy  the  deformation  theory  of  plasticity.  This 

can  be  accomplished  as  follows.  Equation  3»  "the  "k  cycle"  stress 
equation  of  the  preceding  section,  can  be  used  intact  if  it  Is  understood 
(k) 

that  Px  '  is  the  peak  load  at  which  the  results  are  required  and  does  not 
m 

change  from  one  cycle  to  another  as  before.  We  must  iterate  to  a  solution 
in  order  to  obtain  a  satisfactory  approximation  to  the  plastic  strains. 
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The  intra-cycle  procedure  employed  for  the  determination  of  the 

equivalent  strain  for  the  kth  cycle  is  the  same  as  before,  namely  the 
constant  strain  method.  At  this  point,  however,  the  equivalent  strain 
itself,  not  its  increment,  is  resolved  into  the  node  plastic  strains  by 
utilizing  an  engineering  adaptation  of  the  incremental  relations, 
Equations  7>  thus: 


These  are  now  available  for  the  stress  equation  of  the  next  cycle. 

Three  analyses,  one  at  £  =  11,600  lbs.,  one  at  P  =  14,600  lbs., 
and  one  at  P  =  16,760  lbs.  were  performed  on  the  shear  lag  specimen  by 
this  deformation  theory  procedure.  The  results  were  practically  identical 
with  those  shown  in  Figures  (l6)  to  (2l).  The  convergence  to  each  of 
these  results  was  obtained  alter  less  than  thirty  cycles  of  iteration. 

The  machine  time  for  each  calculation  was  approximately  four  minutes. 

In  the  case  of  solutions  like  this,  where  the  two  analyses  give 
practically  identical  results,  the  deformation  approach  is  naturally 
very  attractive  because  of  the  greatly  reduced  machine  time.  However, 
the  question  remains  of  determining  when  to  expect  the  results  to  agree 
in  this  manner. 

The  IBM  pro  gran  presented  in  Appendix  V  cannot  be  used  for  deforma¬ 
tion  theory  analyses.  However,  minor  changes  can  be  made  in  the  program 
to  permit  calculations  of  this  type  to  be  carried  out. 
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SECTION  IV 

ISOTROPIC  ELASTIC-PLASTIC-CREEP  ANALYSIS 


A.  Introduction  to  Creep  Theory 

Strains  due  to  creep  constitute  an  additional  form  of  inelastic  strain, 
and  can  be  handled  in  a  way  analogous  to  that  already  discussed  for  bi¬ 
axial  plasticity  by  the  flow  theory.  It  is  only  necessary  to  select  a 
method  for  evaluating  these  time-dependent  strains  baaed  upon  the  material 
properties  and  to  add  them  to  the  plastic  strains  prior  to  insertion  in 
the  basic  Equation  3* 

A  familiar  relationship  used  to  match  the  creep  behavior  in  a  tensile 
creep  test,  performed  at  constant  stress  and  constant  temperature ,  is 
(Reference  4): 

e(c)  =  atY(ePo-l)  (10) 

in  which 

e(c)  is  the  tensile  creep. strain 

t  is  the  elapsed  time 

a  is  the  constant  tensile  stress 

a,3,Y  are  empirical  constants  for  the  particular  test 
temperature 

For  this  analysis  the  assumption  is  made  that  there  exists  an  effec¬ 
tive  creep  strain  e(c)  in  a  biaxial  situation  which  can  be  calculated 
using  Equation  10.  In  doing  this  the  stress  a  is  taken  to  be  the  von 
Mises  effective  stress  obtained  fran  Equation  5*  The  further  assumption 
is  made  that  this  effective  creep  strain  can  be  resolved  into  nodal  creep 
strains  by  use  of  a  Frandtl-Reuss  type  of  flow  law. 

The  creep  strain  calculation  must  be  generalized  to  situations  in 
which  the  stresses  vary  with  time.  One  well-known  procedure  for  doing 
this,  the  strain-hardening  rule  (Reference  4),  has  been  determined  to  be 
most  appropriate  for  the  present  purposes.  Its  use  will  be  described 

presently  in  connection  with  the  k  calculation  cycle. 

B.  Creep  Theory  Details 

The  calculation  cycle  follows  a  sequence  similar  to  that  described 
previously  for  the  isotropic  elastic-plastic  analysis.  An  additional 
step  is  necessary,  just  after  the  plastic  strains  are  obtained,  to  de¬ 
termine  the  creep  strains.  The  intra-cycle  order  of  calculation  is  a3 
follows. 
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Entering  the  cycle  with  applied  loads  }  and  with  the 

initial  strains,  plastic  and  creep,  as  calculated  during  the  preceding 
cycle,  the  stresses  jo^^)  j  are  calculated  fran  Equation  3» 

The  effective  stresses  obtained  from  Equation  5  ere  used  to 

obtain  the  effective  plastic  strains. 


The  creep  strain  increments  at  each  node,  for  a  specified  time  step, 
are  now  determined  by  the  strain  hardening  rule  which  relates  the  strain 
at  a  node  to  the  corresponding  stress  and  strain  for  the  previous  cycle 
by  the  introduction  of,  an  assumed  elapsed  time. 

Referring  to  Figure  (23),  one  goes  to  the  constant  effective  stress- 

temperature  curve  (ok,  T1*)  relevant  to  the  node  and  the  cycle,  and  locates 

upon  it  the  point  with  ordinate  e(c)(k"^.  The  corresponding  abscissa, 

designated  t*W,  is  called  the  reference  time  and  is  generally  different 
from  the  actual  elapsed  time  at  the  start  of  the  cycle.  The  required 

effective  creep  strain  increment  Ae(c)^  is  that  corresponding  to  the 

increase  in  time  from  t*^  to  (t*^  +  At^)  as  shown  on  Figure  (23),  At^ 
being  the  selected  calculation  time  increment. 


v 

The  increment  A?(c)  is  substituted  into  the  Prandtl-Reuss  type  in¬ 
cremental  relations.  Equations  7>  together  with  the  stresses  indicated 
there.  The  creep  initial  strains  are  then  obtained  as  in  Equation  8. 

In  summary,  the  steps  in  the  k^*1  calculation  take  the  following 
order: 


(1)  Evaluation  of  Equation  3  to  obtaining  the  '•tress  components, 

oOO 

u 

(2)  ®ie  calculation  of  effective  stress  according  to  Equation  5 

(3)  The  determination  of  the  node  plastic  strains 

(4)  The  determination  of  the  node  creep  strains 

(5)  The  addition  of  the  nodal  plastic  and  creep  strains  to  give 
the  initial  strains  for  the  next  cycle. 


C.  Description  of  Structure  and  Tests 

The  description  of  the  shear  lag  structure  to  be  analyzed  in  this 
section  and  tests  for  the ’material  properties  may  be  found  in  Reference  7« 
Hie  shear  lag  structure  was  manufactured  fran  1100-F  aluminum.  It  was 
of  the  same  physical  dimensions  as  the  structure  of  Figure  (6).  The 
Idealisation  of  the  upper  right  quadrant  remains  unchanged. 
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Material  properties  for  the  creep  analysis  were  obtained  from  uni¬ 
axial  strain-time  tests  for  constant  tensile  stress.  The  temperature 
at  which  these  tests  were  conducted  was  206*0,  the  temperature  identical 
to  that  of  the  structural  test.  The  curves  for  these  tests  are  presented 
in  Figure  (24)  together  with  the  fitted  curves  from  Equation  10.  The 
constants  of  the  equation  were  obtained  fran  Reference  7  and  are  as 
follows : 

a  =  0.650  x  10-4 

y  =  0.500 

0  =  0.700  x  10-3  in2/lb 

Ordinary  tensile  stress-strain  tests  were  performed  at  roam  tempera¬ 
ture  on  coupons  cut  from  the  x  and  y  orientations  of  the  plate  material. 
The  data  and  faired  curve  are  presented  in  Figure  (25).  Tensile  stress- 
strain  data  for  206°C  are  also  plotted  here.  These  data  were  obtained 
fran  the  intersections  of  the  test  curves  on  the  zero  time  axis  in  Fig¬ 
ure  (24).  A  piecewise  linear  representation  TCI  was  fitted  to  this 
latter  data. 

The  location  of  the  strain  gages  on  the  structure  is  given  in  Fig¬ 
ure  (26) .  The  shear  lag  specimen  was  tested  for  a  total  of  three  hours 
at  206uC .  An  initial  applied  load  of  1600  lb.  was  increased  to  2020  lb. 
at  the  end  of  the  first  hour.  It  was  held  constant  thereafter  to  the 
end  of  the  test. 

D.  Results  of  Creep  Shear  Lag  Analysis 

The  predicted  strain  distributions  along  the  x-axis  and  along  the 
section  x  =  1  in.  at  t  =  0.06  hr.,  t  =  1.10  hr.,  and  t  =  3. 00  hr. 
elapsed  times,  are  given  in  Figures  (27)  through  (32),  together  with 
the  experimental  data  of  Reference  7«  Test  data  are  not  available  for 
the  y-node  strains  at  the  center  node,  and  so  this  correlation  point  of 
critical  significance  does  not  exist. 

The  curve  of  Figure  (33),  effective  stress  versus  strain  at  the 
central  node,  exhibits  the  shapes  characteristic  to  the  various  regions 
of  the  load-time  sequence.  The  initial  linear  segment ,  representing 
elastic  loading  is  followed  by  the  region  of  the  negative  curvature 
representing  loading  into  the  plastic  range,  all  at  assumed  zero  time. 
Thereafter,  the  applied  load  remains  constant  for  one  hour,  during  which 
time  there  is  a  stress  redistribution  in  the  structure  due  to  creep. 

This  particular  node  unloads,  as  evidenced  by  the  reduction  in  effective 
stress,  although  the  total  strain  is  growing  continuously.  The  applied 
load  is  now  increased  to  2020  lb.  Because  of  the  previous  elastic  strain 
recovery,  the  effective  stress  at  first  goes  up  elastically,  and  then 
becomes  plastic  once  more.  Once  the  applied  load  reaches  its  final  value, 
redistribution  due  to  creep  effects  again  takes  place. 

In  the  initial  stages  of  creep  the  curve  is  very  sensitive  to  time 
increment  size  and  it  is  necessary  to  choose  exceptionally  man  time 


17 


increments  in  this  region  if  good  accuracy  is  to  be  achieved.  The  time 
increments  employed  are  as  shown  in  Figure  (33) • 

Considering  the  simplicity  of  the  expressions  employed  to  describe 
as  complex  a  phenomenon  as  creep  and  the  liberal  assumptions  made  in  the 
process,  the  correlation  between  analysis  and  experiment,  as  evidenced 
by  the  preceding  graphs  and  also  by  Figures  (3*0  and  (35),  is  surprisingly 
good. 
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SECTION  V 


ANISOTROPIC  ELASTIC - PLASTIC  ANALYSIS 


A.  Anisotropy  In  Structure3 

Several  "expanding  yield  surface"  theories  for  extending  the 
isotropic  plastic  theory  to  provide  for  anisotropy  have  became  available 
within  the  last  two  decades.  Each  is  based  on  experimentally  determined 
parameters  and,  therefore,  each  is  biased  in  favor  of  specific  test  data. 
The  complexity  and  amount  of  testing  required  to  obtain  these  parameters 
differ  considerably.  Under  these  circumstances,  no  single  theory  can  be 
completely  acceptable,  but  it  is  thought  that  a  suitable  theory  must,  at 
least,  be  capable  of  evaluating  the  type  of  anisotropy  associated  with 
biaxially  stressed  structures  used  in  flight  vehicle  design  without  being 
unduly  complex  in  application.  It  would  be  desirable  to  have  the  pro¬ 
cedure  based  on  a  well-known,  accepted  theory. 

A  particular  type  of  anisotropy,  the  so-called  "orthotropic  symmetric" 
type,  develops  during  a  cold  rolling  process  where  the  material  iB  length¬ 
ened  and  thinned  with  no  appreciable  change  in  width.  Since  cold  rolled 
sheet  and  plate  are  frequently  used  in  aerospace  structures,  this  type  of 
anisotropy  may  be  anticipated  and  is  considered  here. 

A  theory  proposed  by  Hill,  Reference  10,  has  been  widely  accepted 
as  the  most  straightforward  extension  of  the  isotropic  theory.  The  formu¬ 
lation,  however,  is  not  very  convenient  for  numerical  step-by-step  compu¬ 
tation.  A  modification  of  Hill's  theory  proposed  by  Hu,  Reference  11, 
however,  is  very  tractable  to  formulation  into  the  matrix  inelastic  pro¬ 
gram  discussed  previously  in  this  report.  The  Hu  procedure  has  two  dis¬ 
tinct  advantages : 

(1)  It  employs  a  von  Mises  type  hardening  surface,  associated  flow 
law  and  effective  stress-strain  relationships  in  appropriate 
foim. 

(2)  It  requires  a  minimum  of  material  data:  simple  uniaxial  and 
shear  stress-strain  tests  on  coupons  cut  in  the  directions  of 
the  orthotropy. 

B.  Hu's  Strain  Hardening  Theory 

A  summary  of  Hu's  theory  is  presented  to  establish  its  limitations 
and  provide  background  for  the  necessary  modifications  to  obtain  a  more 
general  theory. 

The  isotropic  expressions  for  effective  stress,  Equation  5,  and 
associated  incremental  relations,  Equations  7,  are  modified  by  means  of 
anisotropic  parameters  (or^).  These  are  constants  in  Hu's  theory.  Here 

it  is  more  convenient  to  introduce  the  1-1,  2-2,  etc ., directions  Instead 


of  x-x,  y-y,  etc.  The  modified  expression  for  the  effective  stress  is 
o  =  [a^  an  -  a22 )  +  o>23(a22  -  o^)  +  «31(  a33  -  CTn) 

L  i  (11) 

+  +  3or55°23  +  ^e^h]2 

The  incremental  flow  equations  are 
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Equations  11  and  1 2  are  written  for  the  case  where  the  reference  axes 
are  the  principal  axes  of  anisotropy. 

The  anisotropic  parameters  are  determined  by  means  of  a  total  of 
six,  simple,  directional,  stress-strain  tests  (l.e.,  uniaxial  and  shear 
tests),  where,  alternately,  all  stress  components  are  equal  to  zero 
except  one.  Ercm  each  of  the  six  tests  a  characteristic  stress,  such 
as  an  approximate  yield  stress,  is  read  off.  Then  substituting  each 
of  these  results  into  Equation  11  in  succession,  we  may  write 


an  =  aL2  +  “31  “  (^1) 

“22  =  “23  +  ®L2  =  (^) 

»55-tfe)a 

a33  =  “31  +  a23  =  (^) 

where 


K  =  the  effective  characteristic  stress 

Ty  *  the  simple  directional  characteristic  stress 


It  remains  to  assume  an  effective  stress-effective  strain  relation¬ 
ship.  Hu  shows  that  it  is  acceptable  to  assume  the  stress-strain  data 
associated  with  one  particular  simple  tension  test  as  the  effective  re¬ 
lationship.  It  can  be  seen  fran  Equations  13  that  this  implies  setting 
one  of^  equal  to  unity. 


Thus,  in  particular  for 

=  1,  from  Equations  11 

end  12: 

5 

de 
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W 

Now  consider  Equation  12  for  simple  stress-strain 
directions.  Then 
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These  equations  3ay  that,  with  the  anisotropic  parameters  constant  for 
strain  hardening,  the  simple,  stress-incremental  strain  relations  must 
be  proportional  to  the  effective  stress-incremental  strain  relationship. 
The  implication  is  that  the  integrated  forms  of  Equations  15,  that  is, 
the  simple,  directional  stress  strain  curves,  are  thus  prescribed.  These 
may  or  may  not  be  a  reasonable  fit  to  the  test  data  for  the  material  Of 
interest.  Obviously,  only  when  the  fit  is  good  can  one  hope  for  Hu's 
theory  to  give  acceptable  results  for  all  types  of  loading. 

Based  upon  the  Hu  theory,  it  becomes  relatively  easy  to  obtain 
anisotropic  solutions  using  the  previously  developed  isotropic  inelastic 
procedure  and  corresponding  digital  computer  program.  It  is  only  nec¬ 
essary  to  substitute  the  appropriate  anisotropic  constants  for  their 
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isotropic  counterparts,  for  which 


**“66  "  x>  *12  =  *23  =  “31 


1 

2* 


C.  Extension  of  Hu's  Theory 


It  is  in  the  determination  of  the  parameters  that  an 'extension  to 
Hu's  theory  is  made.  Except  for  the  special  case  pointed  out  in  the  pre¬ 
ceding  section,  these  parameters  should  not  be  constant  for  a  strain 
hardening  material ^  but  should  be  variables  dependent  upon  stress  level.' 
The  objective,  obviously,  is  to  determine  the  variation  in  a  manner  that 
allows  for  all  of  the  simple,  directional  stress-strain  curves  to  be 
correctly  reproduced.  This  can  be  accomplished  by  a  consideration  of 
plastic  work.  .  ,  . 

In  the  current  approach,  one  continues  with  the  assumption  of  the 
existence  of  an  effective  stress-strain  relationship.  Then  the  basic 
notion  is  that  the  anisotropic  parameters  are  determined  such  that,  for 
equal  amounts  of  plastic  work  done  during  simple  directional  stress 
strain  tests  in  all  directions,  the  effective  stress  level  reached  will 
be  identical. 

I 

Accordingly,  for  a  tensile  specimen  in  the  1-1  direction,  one  calcu¬ 
lates  the  plastib  work  w  performed  during  a  uniaxial  test  by  the  formula 


w 


Ai  den  =  A  ** 


(16) 


Let  the  corresponding  maximum  stress  reached  be  identified  by  the  super¬ 
script  (I),  i.e.  <£.  For  a  similar  test  in  the  2-2  direction,  and  for 
which  the  amount  of  plastic  work  performed  is  identical,  the  corresponding 


maximum  stress  reached  is  a. 
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Since  the  amounts  of  work  done  in  the  two 


I  I 

cases  are  the  same,  a  and  a  correspond  to  the  same  effective  stress 

j  XI 

a.  By  means  of  Equation  11  we  have 


(/!  ■  ■  v,(4)S 


(17) 


This  expression  constitutes  a  relationship  defining  7  and  as 

functions  of  5.  Similar  relationships  can  clearly  be  found  for  the  ether 
c^'s.  Thereafter,  the  c^'s  can  be  determined  as  functions  of  o  by  re¬ 
course  to  the  a  definitions  of  Equations  13 . 

It  is  convenient  to  again  select  the  2-2  direction  stress-strain 
curve  as  the  one  defining  the  effective  stress-strain  relationship. 

This  results  in  being  equal  to  unity  once  more. 
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The  actual  evaluation  of  the  other  a's  as  functions  of  3  follows 
easily.  Figure  (3 6)  indicates  schematically,  how  the  plastic  work  done 
in  each  simple  directional  stress-strain  test  can  be  plotted  as  a  function 
of  stress.  Then  for  a  given  amount  of  plastic  work,  and  reading  off  the 

corresponding  stresses  c*  and  o*  for  example,  by  Equation  IT  one  finds 
that 


Using  this  and  similai*  relationships  in  the  other  directions,  curves 
.representing  all  of  the  a's  as  functions  of  .8  may  be  constructed.  The 
incorporation  of  this  information  in  the  step-by-step  calculation  pro¬ 
cedure  is  discussed  in  the  next  section. 


The  detailed  step-by-step  calculation  procedure  to  be  followed  in 
the  case  of  anisotropic  material  in  a  biaxially  stressed  structure  is 
very  similar  to  that  previously  discussed  for  isotropic  materials. 
Accordingly,  only  the  differences  will  be  stressed. 

As  before,  one  starts  the  kth  calculation  cycle  by  evaluating 
[o^ }  by  means  of  Equation  3  •  This  operation  employs  the  initial 
strains  of  the  preceding  cycle 


The  next  operation  is  to  evaluate  the  effective  stresses  at  each 
of  the  nodes,  3^  ,  by  Equation  11,  modified  for  the  biaxial  case  to 


We  continue  here  our  assumption  that  the  2-2  direction  has  been  selected 
as  that  in  which  the  effective  stress-strain  relationship  is  defined; 
hence  a^  *  !• 
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Note  that  if  the  Hu  theory  ie  being  used,  the  o^'s  are  all  known 

constants.  If  the  modified  theory  is  being  employed,  the  variable  or^'s 

are  those  that  have  been  evaluated  during  the  preceding  cycle,  (k-l). 
This  is  in  keeping  with  the  overall  nature  of  the  analysis  as  a  step-by- 
step  procedure . 

(k) 

Having  determined  the  '  s,  one  can  now  go  to  the  curves  repre- 


senting  the  or^'s  as  functions  of  5  to  evaluate 


00  (k) 


12 


and  orj^ . 


After  having  also  determined  Sjj(p)^  and  AijjCp)^  as  before  by  the 

constant  strain  method,  these  quantities  msy  be  substituted  into  a 
finite  equivalent  of  Equations  12  specialized  for  the  biaxial  case,  to 
yield 
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Fran  these,  the  strain  components  c(p)^k^,  c(p) and  e(p)^00  are 
obtained  by  addition,  as  before,  using  Equation  8. 


After  incrementing  the  applied  load  the  sequence  can  now  be  re¬ 
peated  for  the  next  load  cycle. 


E.  Rotation  of  Axes  of  Anisotropy 

Anisotropic  symmetry  may  occur  in  a  structure  for  which  it  is  con¬ 
venient  to  choose  coordinate  axes  that  are  rotated  from  the  orthogonal 
axes  of  anisotropy.  The  corresponding  expressions  for  the  incvenentnl 
flow  equations  and  effective  stress  equations  are  derived  in  Appendix  17. 
The  derivation  is  limited  to  the  case  of  biaxial  stress  where  the  3-3 
and  z-z  axes  coincide.  A  method  of  obtaining  the  shear  anisotropic  co¬ 
efficients  is  also  indicated. 
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P.  Anisotropic  Analysis  of  Shear  Lag  Structure 

It  has  been  pointed  out  earlier  that  the  202k-lb  aluminum  alloy  of 
the  shear  lag  structure  tested  for  the  Air  Force  (Reference  7)  displayed 
considerable  anisotropy  as  shown  in  Figure  (ll).  The  structure  had  "been 
analyzed  for  isotropic  strain  hardening  based  upon  the  curve  R034  in 
this  figure,  and  the  results  discussed  in  Section  III,  A  corresponding 
anisotropic  strain  hardening  analysis  has  also  been  carried  out,  employing 
first  the  Hu  theory  and  then  the  proposed  extension. 

Material  stress-strain  data  is  available  along  two  of  the  axes  of 
anisotrow,  the  rolling  or  x-x  (l-l)  direction,  and  the  long  transverse  or 
y-y  (2-2)  direction.  This  data  has  been  plotted  and  discussed  in 
connection  with  Figure  (ll);  it  is  replotted  for  convenience  in  Figure  (37)* 

The  additional  required,  but  missing,  test  data  is  (a)  tensile  stress- 
strain  data  in  the  short  transverse  or  z-z  (3-3)  direction  and  (b)  shear 
stress-strain  data  in  the  x-y  plane.  A  reasonable  assumption  to  make  for 
engineering  purposes  for  (a)  is  that  the  long  and  short  transverse  properties 
are  identical;  this  is  made  in  the  analyses  to  follow.  In  the  case  of  (b), 
the  missing  shear  data,  the  following  is  done.  First,  a  shear  stress-strain 
curve  is  obtained  based  upon  the  tensile  curve  in  the  rolling  direction, 
together  with  the  assumption  that  the  material  is  isotropic  and  governed  by 
the  incremental  theory  of  Section  III  —  specifically.  Equations  5  and  7. 
Next,  a  similar  shear  stress-strain  curve  based  upon  the  tensile  curve  in  the 
long  transverse  direction  is  obtained.  Finally,  a  faired- in  average  of  these 
two  curves  is  taken  to  represent  the  missing  shear  stress-strain  relation. 

In  order  to  apply  the  J*a  theory,  one  must  first  select  four  character¬ 
istic  stresses  to  represent  the  directional  stress-strain  curves,  as  dis¬ 
cussed  in  Section  V-B.  These  are  the  quantities  Yg2,  Y^  and  Y12  of 

Equations  13 .  We  arbitrarily  choose  the  proportional  limits  from  the  four 
curves  Just  discussed  for  these  values;  they  are 

Yn  =  51  ksi 
Y22  -  32  ksi 
Y33  =  51  ksi 
=  22  ksi 


Since  the  simple  directional  stress-strain  curve  in  the  2-2  direc¬ 
tion,  i.e.  R02M,  is  selected  for  the  effective  stress  relationship,  the 
effective  stress  characteristic  value  K  is  also  32  ksi.  Substituting 
these  values  in  Equations  13,  one  can  solve  for  the  necessary  anisotropic 
parameters.  This  infoimation  is  nil  that  is  required  to  carry  out  an 
analysis  of  the  shear  lag  structure  by  the  Hu  theory. 

As  pointed  out  previously,  once  the  anisotropic  parameters  cr.  ,  a 

1  1.  C.C.  j 

o/j2  and  Qf}|  t)  have  been  specified  and  one  of  the  stress-strain  curves 
chosen  as  'the  effective  stress-strain  curve,  the  remaining  simple  direc- 
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tional  stress-strain  curves  are  prescribed.  Let  us  new  examine  the  con¬ 
sequences  of  our  current  parameter  selection.  Since  only  test  data  in 
the  1-1  and  2-2  directions  are  available,  ve  shall  concentrate  on  these. 
For  this  special  case,  Equations  11  and  12  can  be  easily  manipulated  to 
yield 


5  =  /gn°ii 


den 

do-Li 


de 

an  dd 


(21) 


Using  these  expressions,  one  can  construct  a  stress-strain  curve. 

Such  a  curve  is  shown  on  Figure  (37)*  As  can.be  seen,  the  fit  with 
the  1-1  test  data  is  very  poor.  Apparently,  one  could  do  better  by 
arbitrarily  choosing  for  the  1-1  characteristic  value  a  much  lower 

value  than  51  kai  —  perhaps  in  the  neighborhood  of  35  ksi.  Nevertheless, 
it  is  clear  that  in  this  case  it  would  still  be  impossible  to  get  a 
really  good  fit  because  of  the  fundamentally  different  shapes  of  the  two 
stress-strain  curves,  especially  in  the  region  of  their  knees.  It  will 
turn  out  that,  in  the  case  of  the  shear  lag  structure,  this  selection 
is  not  critical,  because  of  the  fact  that  stresses  in  the  1-1  direction 
are  very  low,  compared  to  those  in  the  2-2  direction. 

In  order  to  carry  out  an  analysis  of  the  shear  lag  structure  based 
upon  the  extension  to  the  Hu  theory,  one  must  first  evaluate  the  plastic 
work  done  in  each  of  the  simple  directional  stresss-strain  tests  as  a 
function  of  the  applicable  stress.  Using  this  information  in  the  manner 
discussed  previously,  one  can  then  obtain  the  anisotropic  parameters  as 
functions  of  effective  stress.  This  has  been  done,  assuming  that  it  is 
sufficiently  accurate  to  represent  the  curves  by  a  small  number  of 
connected  straight  line  segments.  The  key  values  of  the  resulting  cr’s 
are  given  in  the  accompanying  table.  Corresponding  total  strains  £ 

(  =  a/E  +«)  are  also  listed. 
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In  the  preceding  table,  the  reduction  number  of  independent  anisotropic 
parameters  is  a  consequence  of  the  assumptions  which  had  to  be  made  for 
the  missing  stress-strain  curves. 


The  value  of  unity  for  *n  &ccordance''vith  our  continued 

selection  of  the  tensile  stress-strain  curve  R02M  incite  y-y  or  2-2 
direction  as  the  basis  for  our  effective  stress-strain  relationship. 

It  might  be  reiterated  that,  because  of  the  manner  in  vhitet^the  or’s 
are  derived  in  this  case } all  of  the  simple  directional  stress- strain 
curves  are  matched  as  closely  as  desired. 


Using  first  the  anisotropic  parameters  selected  as  described  for 
the  Hu  theory,  and  then  for  its  extension,  the  shear  leg  structure  has 
been  analyzed.  The  results  are  presented  in  the  next  section. 


G.  Discussion  of  Results 

Sane  of  the  highlights  of  the  two  anisotropic  analyses  are  pre¬ 
sented  in  Figures  (38)  through  (4o) .  Corresponding  results  obtained 
previously  and  based  upon  isotropic  theory  are  also  included  for  com¬ 
parison. 

Figure  (38)  refers  to  the  central  node  of  the  shear  lag  structure 
and  shows  the  variation  of  the  total  strains  in  the  x  and  y  directions 
with  increasing  applied  load.  The  isotropic  result  is  replotted  from 
Figure  (22)  specifically  those  curves  based  upon  the  effective  stress- 
strain  curve  R031.  ' 

It  can  be  seen  from  Figure  (38)  that  the  most  flexible  analysis 
predictions,  that  is,  those  for  which  the  total  strains  are  the  largest, 
are  obtained  by  the  isotropic  analysis.  The  analysis  based  upon  the 
extension-  to  the  Hu  theory  is  somewhat  less  flexible,  while  the  Hu 
theory  results  are  the  stiffest.  Also,  the  differences  between  the 
three  analyses  are  less  in  the  y  direction  than  in^the  x  direction. 

These  observations  are  clearly  in  agreement  with  the  nature  of 
the  material  stress-strain  relations  upon  which  the  analyses  are  based. 

In  alt  three  cases,  we  are  assuming  that  the  y-y  curves  are  identical, 
that  is,  all  are  represented  by  R02M.  Referring  to  Figure  ( 37) ,  in  the 
case  of  the  isotropic  analysis,  we  are  also  assuming,  in  effect,  that 
the  x-x  stress-strain  curve  is  identical  to  the  y-y  curve  R02M.  The 
x-x  curve  for  the  extension  to  the  Hu  theory  is  seen  to  be  somewhat 
stiffer  than  R02M,  while  the  x-x  curve  for  the  Hu  theory  is  by  far  the 
stiffest  of  the  three. 

Comparing  the  analyses  with  the  test  data,  all  three  analyses  sub¬ 
stantially  underpredict  the  test  points  in  the  high  applied  load  regime, 
and  more  so  in  the  y  than  in  the  x  direction.  It  is  interesting  to  note 
however,  that  while  the  anisotropic  analyses  make  these  differences  even 
greater,  they  do  have  the  virtue  of  making  the  comparison  more  consistent 
as  between  the  x  and  y  directions.  Thus  it  would  appear  that  the 
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anisotropic  nature  of  the  2024-T4  material  does  influence  the  strain 
distribution  in  the  shear  lag  structure,  and  that  the  anisotropic  analyses 
can  detect  this  tendency. 

Figures  (39)  and  (40)  show  the  strain  distributions  along  the  two 
strain  gage  lines  of  the  test  structure.  In  addition  to  the  isotropic 
predicted  results,  replotted  frcm  Figures  (20)  and  (2l),  anisotropic 
results  based  upon  the  extension  to  the  Hu  theory  are  presented.  They 
indicate  that  while  the  differences  are  not  dramatic,  they  do  in  fact 
exist . 
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SECTION  VI 


ANISOTROPIC  ELASTIC-PLASTIC-CREEP  ANALYSIS 


A.  Anisotropic  Creep  Theory 

The  application  of  matrix  analysis  procedures  to  problems  involving 
anisotropic  creep  is  at  present  academic.  Neither  meaningful  analytical 
research  nor  appropriate  test  data  has  been  found.  Therefore  the  pro¬ 
cedure  presented  here  for  anisotropic  creep  is  a  simple  extension  of 
that  already  described  in  previous  sections.  Accordingly,  only  the 
difference  and  additional  assumptions  are  discussed. 

Strains  due  to  anisotropic  creep  can  be  handled  in  a  way  similar  to 
those  due  to  isotropic  creep.  The  step-by-step  procedure  has  only  to  be 
modified  for  the  anisotropic  behavior  by  substituting  Equations  19  and 
20  for  Equations  5  and  T  respectively  in  the  manner  described  for  time 
independent  plastic  anisotropy  in  Section  V-D.  The  additional  assumptions 
implied  by  this  simple  extension  are  the  following: 

(1)  The  anisotropic  parameters  calculated  fran  zero  time  simple 
directional  tests  (by  either  Hu  theory  or  the  proposed  extens¬ 
ion)  are  valid  for  anisotropic  creep. 

(2)  The  effective  creep  strain  equation,  Equation' 10,  remains  valid 
and  the  empirical  constants  (a,fj,y)  are  determined  for  the 
tensile  creep  test  in  the  assumed  effective  stress-strain  di¬ 
rection. 

Because  testing  of  structures  exhibiting  anisotropic  creep  has 
apparently  not  been  done,  the  anisotropic  creep  procedure  cannot  be 
checked  out  against  teats  at  this  time.  However,  the  results  of  a 
sample  calculation  are  presented  for  a  hypothetical  material  having 
this  characteristic. 

B.  Sample  Problem 

The  1100-F  aluminum  shear  lag  specimen,  already  analyzed  for  isotropic 
creep,  is  used  for  an  anisotropic  creep  analysis  making  the  following 
assumptions: 

(1)  All  the  uniaxial  data  (Figure  25)  employed  in  the  isotropic 
creep  analysis  is  assumed  to  refer  to  the  y-y  (2-2)  effective 
direction  of  the  anisotropic  creep  analysis. 

(2)  The  plastic  anisotropic  parameters  (extension  of  Hu's  theory) 
for  the  effective  stress-strain  curve  of  the  2024-T4  aluminum 
alloy  material  for  successive  levels  of  stress,  32  ksi,  37  ksi, 

40  ksi,  etc.,  are  arbitrarily  chosen  for  the  stress  levels, 

2  ksi,  3  ksi,  4  ksi,  etc., of  the  1100-F  aluminum  effective 
stress-strain  curve  TCI  (Figure  25).  Thus: 
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This  means  that  the  uniaxial  x-x  (l-l)  predicted  stress-strain  curve  and 
the  shear  x-y  stress-strain  curve  are  arbitrarily  stiffer  than  the 
corresponding  isotropic  curves. 

•  The  node  strains  obtained,  using  the  anisotropic  assumptions  for 
load  P  =  2020  lbs  at  time  t  =  3  hr s,  are  presented  for  the  x-axis  and 
along  x  =  1  in.  in  Figures  (4l)  and  (42)  respectively. 

In  general,  the  anisotropic  results  appear  to  be  stiffer  than  those 
for  the  isotropic  case,  and  this  is  especially  true  for  the  vertical  gage 
line,  where  the  anisotropic  shear  strains  are  of  the  order  of  3/4  of  the 
isotropic  values.  This  should  be  expected,  due  to  the  fact  that  the 
anisotropic  parameters  used  are  all  less  than,  or  at  the  most  equal  to, 
their  isotropic  equivalents.  '  ' 
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CONCLUSION 


The  linear,  matrix  structural  analysis  methods  currently  in  general 
use  throughout  the  aerospace  industry  have  recently  been  extended  to  in¬ 
clude  structures  loaded  into  the  inelastic  material  behavior  regime.  However, 
very  little  published  information  is  available  correlating  predicted 
results  with  the  test  data. 

The  current  report  recanmends  a  simple  analytical  approach  to  such 
problems.  It  is  based  upon  the  concept  of  initial  strains  in  combination 
with  a  suitable  matrix  of  influence  coefficients,  obtained  by  standard 
linear  matrix  structural  analysis  methods.  The  initial  strains  are  those 
associated  with  plasticity  and  creep. 

Sane  particularly  useful  tests  of  aluminum  and  aluminum  alloy  shear 
lag  structures  have  been  performed  previously  for  the  Air  Force.  These 
structures  have  been  analyzed  by  the  recommended  method,  and  the  result¬ 
ing  agreement  (for  both  the  plasticity  and  the  creep  tests)  is  considered 
to  be  very  encouraging.  On  the  other  hand,  additional  testing  must  be 
carried  out  and  correlations  with  analysis  made  before  the  method  can 
be  considered  as  fully  evaluated. 

In  the  meantime,  the  writers  believe  that  sufficient  confidence  in 
the  method  has  been  established  that  it  may  now  be  used  in  practical 
engineering  applications.  For  example,  it  should  be  immediately  useful 
in  such  problems  as  predicting  inelastic  strain  distributions  around 
stress  raisers  in  simple  structural  canponents. 

As  for  the  influence  of  anisotropy,  it  has  been  shown  that  this 
material  property  can  be  readily  accommodated  in  the  recommended  procedure. 
Calculations  made  for  the_  2024-Tit-  shear  lag  structure  indicate  that,  in 
this  case  at  least,  anisotropy  plays  a  minor  but  discernable  role  in  de¬ 
termining  the  strain  distributions.  This  evidence  is  inconclusive, 
because  in  this  particular  test,  while  the  material  itself  is  decidely 
anisotropic,  the  stresses  normal  to  the  direction  of  the  applied  loads 
are  quite  small. 

It  is  recommended  that  in  cases  where  doubt  exists  as  to  the  im¬ 
portance  of  anisotropy,  it  be  included  in  the  analysis.  Certainly  in 
such  examples  as  the  shear  lag  structure,  the  additional  complexity  in 
the  use  of  the  computer  program  is  very  slight. 

There  is  one  very  important  restriction  implicit  in  the  method  pro¬ 
posed  in  this  report.  Essentially,  this  method  applies  only  to  a  struc¬ 
ture  in  which  the  material  is  initially  in  the  virgin  state,  and  there¬ 
after  experiences  only  continually- increasing  applied  loads  until  failure 
occurs .  This  limitation  can  actually  be  relaxed  to  the  extent  that 
elastic  unloading  followed  by  reloading  can  also  be  accommodated,  but 
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completely  reversed  loading  is  specifically  excluded.  For  the  latter 
case,  a  different  plasticity  theory  is  needed.  It  is  important  that  a 
technique  for  handling  such  problems  he  developed  because  of  the  need 
for  such  applications  in  fatigue  work.  Efforts  toward  removing  this  re¬ 
striction  are  currently  under  way  under  government  contract. 

Other  than  the  preceding,  the  largest  remaining  obstacle  to  com¬ 
plete  inelastic  analysis  of  practical  aerospace  structures  is  believed 
to  be  the  dearth  of  appropriate  material  property  data,  and  constitutive 
laws  to  describe  them.  1 
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Stress 


Total  Strain 


Fig.  2  Constant  Stress  Method 


Total  Strain 

Fig.  3  Constant  Strain  Method 


Fig.  4  Results  of  Constant  Stress  Method  for  Bar  No.  3  of  Truss 


Load  -  kips 

Fig.  5  Results  of  Constant  Strain  Method  for  Bar  No.  3  of  Truss 
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t  (Thickness) 
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y  -  in. 


Fig.  15  Elastic  Strains  Along  x  =  0.8125  in.  for  P  =  1  lb 
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Fig.  15  Strains  Along  x  ■  0.8125  in.  for  P  ■ 


4 

14,600  lb  and 
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Strain  in/in  x  10-'5  Strain  -  in/ln  x  10 


Fig.  21  Strains  Along  x  =  0.8125  in.  for  P  =  16,7^0  lb  and  AP  *»  5  lb 


Fig.  24  1100-F  Aluminum  Time-Lependent  Behavior  and  Fitted  Curves  (Ref.  7) 


Fig.  27  Strains  ..long  x-Axis  for  P  =  1600  lb  and  t  =  0.06  hr 


Fig.  28  Strains  Along  x  =  1  in.  for  P  =  1600  lb  and  t  =  0.06  hr 


Fig.  32  Strains  Along  x  ■  1  in«  for  P  »  2020  lb  and  t  »  3»0  hr 


Fig.  32  Strains  Along  x  *  1  in*  for  P  •  2020  It  and  t  »  3*0  hr 


Effective  Stress  vs  Strain  at  Center  Mode 


Strains  -  in. /in 


10 


Time  -  hr 


Fig.  34  Total  (Elastic,  Plastic  and  Creep)  Strains 

at  Center  of  Specimen  (Gage  l)  for  P*l600  lb 
to  Time  1  hr,  then  P=2020  lb  to  Time  3  hr 
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Uniaxial  Stress  -  Plastic 
Strain  Showing  Plastic 
Work  w*  =  /  or  dc 
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Tensile  Strain  -  in. /in.  x  10 


Fig.  37  Assumed  Effective  Stress  (y-y)  Curve, 
Calculated  x-x  Curves  and  feat  Data 
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Fig.  39  Strains  Along  x  -  Axis  for  P  =  16,760  lb. 
Isotropic  and  Anisotropic  Analyses 
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Hu's  Theory 
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Fig.  40  Strains  Along  x  =  0.8125  in.  for  P  =  16,760  lb 
Isotropic  and  Anisotropic  Analyses 


Strain  -  in/in  x  10 


x  -  in. 


Fig.  4l  Strains  Along  x  -  Axis  for  P  =  2020  lb  and  t  =  3*0  hr. 
Isotropic  and  Anisotropic  Analyses 


y  -  in. 

Fig*  42  Strains  Along  x  =  1  In.  for  P  =  2020  lb  and  t  =  3*® 
Isotropic  and  Anisotropic  Analyses 
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Fig.  ^3  Elastic  Stress  Distribution  Along  x-Axis 


APPENDIX  I 


ELASTIC  ANALYSIS  OF  SHEAR  LAG  STRUCTURE 


A.  Idealization  of  Shear  Lag  Structure 

As  stated  in  the  introduction,  when  a  problem  is  formulated  by 
means  of  a  standard  influence  coefficient  approach,  the  necessary  linear 
analysis  may  be  carried  out  using  either  the  force  or  displacement 
method.  Since  published  correlations  between  results  of  the  matrix  force 
and  direct  stiffness  methods  of  linear  elastic  analyses  for  redundant 
structures  have,  in  the  past,  left  room  for  doubt  as  to  the  equivalence 
of  results,  this  Appendix  presents  a  comparison  of  the  stresses  frcm  two 
idealizations  of  the  simple  shear-lag  stiff ened-plate  structure,  Figures 
(7)  and  (8).  In  the  past,  discrepancies  have  been  due  in  part  to  a 
marked  difference  in  the  arrangement  of  node  points  for  corresponding 
idealizations,  and  also  to  the  fact  that  techniques  for  obtaining  node 
stresses  in  finite  element  analyses  are  still  being  improved.  An  attempt 
was  made  to  keep  the  idealizations  as  comparable  as  possible  with  respect 
to  location  of  nodes  and  the  determination  of  stresses. 

B.  Force  Method 

The  idealization  for  the  force  method  may  be  seen  in  Figure  (8). 

It  comprises  conventional  bars  and  shear  panels  located  in  the 
manner  shown.  The  analyses  of  sane  previous  idealizations  of  this  type 
have  omitted  the  Poisson's  ratio  effect.  This  effect  can  be  incorporated 
in  the  manner  described  in  Appendix  II. 

C.  Stiffness  Method 

The  idealization  for  the  stiffness  method  consists  of  "Turner 
triangles,"  which  are  located  as  shown  in  Figure  (7).  The  basic  theory 
of  the  triangle  is  to  be  found  in  References  8  and  9- 

While  the  conventional  procedure  was  used  to  obtain  node  stresses 
for  the  force  method,  comparable  stiffness-method  stresses  can  be  cal¬ 
culated  in  several  ways.  A  recent  paper,  Reference  9>  suggests  two  means 
of  obtaining  node  stress,  one  of  which  was  employed  in  the  analysis. 

This  method  will  be  reviewed  briefly. 

Figure  (10)  represents  a  cluster  of  triangles.  It  is  required  to 
find  the  stresses  at  the  node  common  to  triangles  P  to  W.  The  node 
forces  for  each  triangle  at  this  apex  are  obtained  as  described  in 
Turner's  former  papers.  Summing  the  forces  on  a  vertical  section  through 
1  in  both  directions  gives 


X_ 


XP+XQ  + 


YP  +  YQ  +  YR  +  YS 
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Analogous  results  X,  and  Y,  are  obtained  for  the  corresponding  horizontal 
cut.  *h  h 


The  node  stresses  at  the  node  are 


4  (ct  +  dt  ) 
*  p  s 


y  "  4  (at  +  bt 


xy 


-  J  5i  .  Yv  1 

L*  (atr  .  btj  i  (ct  ♦  dt  )J 


where  a,  b,  c,  and  d  are  as  indicated  and  t  is  the  thickness. 


D.  Comparison  of  Elastic  Results  and  Perspective 

The  results  are  correlated  by  means  of  the  curves  appearing  in 
Figures  (43)  and  (44)  with  calculated  values  from  the  experimental  data 
of  Reference  J. 

The  correlation  between  the  stresses  derived  from  the  force  method 
and  from  the  stiffness  method  of  analysis  is  excellent  and  may  be  re¬ 
garded  as  exact  for  engineering  applications.  The  largest  discrepancy 
is  in  the  direct  stresses  in  the  x  direction  at  the  middle  of  the  plate 
as  shown  in  Figure  (43).  Even  in  this  region  the  difference  is  quite 
small.  It  is  believed  that  an  even  closer  agreement  could  be  obtained 
by  modifying  the  idealized  structure  to  provide  square  shear  panels 
adjacent  to  the  reinforcement  and  a  finer  grid  at  the  plate  center. 

The  largest  discrepancy  between  analysis  and  test  results  is 
located  in  the  region  of  the  plate  center.  Reference  7  indicated  that 
considerable  bending  was  exhibited  by  the  structure  as  the  ends  of  the 
stiffener  were  loaded.  The  extent  to  which  this  affects  the  gage 
readings  was  not  determined;  however,  it  may  be  anticipated  that  the 
effect  be  greatest  near  the  middle  of  the  plate.  The  curves  on  Figure 
(43)  reinforce  this  impression. 

On  the  basis  of  the  excellent  agreement  noted  here,  it  can  be  con¬ 
cluded  that  the  incorporation  of  plastic  and  creep  effects  into  the 
present  method  of  structural  analysis  will  not  be  restricted  in  any 
way  by  the  particular  linear  analysis  method  employed. 
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APPENDIX  II 


POISSON'S  RATIO  EFFECT 

i 


The  strain  energy  relationship  for  an  elastic  plate  in  terms  of  in 
plane  stresses  is  given  by  the  volume  integral: 


i  fb 


-  2v  °xxayy  +  +  v^axy] 


The  idealized  structure  corresponding  to  a  rectangular  plate  for 
a  finite  element  matrix  force  analysis  has  orthogonal  bars  taking  only 
normal  stresses  and  panels  in  pure  shear*  Figure  (8),  the  shear-lag 
specimen,  represents  a  typical  idealization  of  this  type.  Defining 


a  and  a  as  the  axial  stress  in  the  bars  and  o 
xx  yy  xy 


as  the  shear  stress 


in  the  panels,  the  strain  energy  U'  of  the  idealized  structure  is 


sometimes  taken  as: 


Comparing  the  plate  strain  energy  U  with  the  idealized  structure 

strain  energy  U'  it  is  obvious  that  the  finite  element  expression  ne¬ 
glects  to  account  for  the  Interaction  term  of  the  plate  (-  2v  °xx°yy) 

which  is  due  to  Poisson's  ratio.  This  uncoupling  of  the  normal  stresses 
has  the  effect  of  making  the  idealized  plate  less  rigid  than  the  actual 
plate.  The  finite  element  idealization  is  refined  by  including  the  term 

(-  2v  a  5  )  in  U'  making  the  model  more  consistent  with  the  elastic 

plate.  **yy  > 

The  sketch,  Figure  ($?)>  shows  a  shear  panel  with  adjacent  axial 
load  carrying  bars.  Assume  the  structure  represents  a  portion  of  a 
plate  structure  /  inches  long,  b  inches  wide  and  t  inches  thick.  The 
normal  stresses  at  one  corner  of  the  idealized  structure  are  designated 
and  Og  in  the  x  and  y  directions  respectively.  It  is  sufficiently 

accurate  in  accounting  for  the  interaction  term  to  assume  that  the  normal 
stresses  are  constant  over  the  plate  corresponding  to  the  shaded  quadrant 
and  also  to  assume  that  these  stresses  are  equal  to  and  Og,  the  values 

at  the  corner.  The  strain  energy  term  to  be  included  is  represented  by 


i  ala2  2  a12 
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where  represents  the  interaction  flexibility  influence  coefficient. 

Carrying  out  the  integration  over  the  quadrant  of  the  plate,  the  in¬ 
fluence  coefficient  is  evaluated. 


v  b/t 

"IT 


This  term  together  with  the  reciprocal  term  and  similar  terms  for 

other  biaxially  stressed  areas,  when  included  in  the  flexibility  matrix, 
account  for  the  Poisson's  ratio  effect. 


APPENDIX  IH 


INSTABILITY  ANALYSIS  OF  THE  CONSTANT  STRESS  AND  CONSTANT  STRAIN  METHODS 


The  stability  of  the  two  methods  is  easily  tested  by  determining 
whether  perturbations  introduced  into  the  analysis  grow  or  decay  with 
succeeding  steps.  It  is  instructive,  however,  to  take  the  following 
approach.  In  place  of  the  computed  quantities,  introduce  exact  quanti¬ 
ties,  signified  by  a  caret  and  associated  error  terms,  in  the  manner 


Wk)M^Me,s- 


III-l 


The  exact  relation  is  taken  to  be 


hi  ■  MKi  *  y  hi  • 


m-2 


Equation  3  can  be  written  in  terms  of  member  loads 

{4kl}  -  *  [qjf'*-1’} 

Substituting  Equation  III-l  into  Equation  III-3  gives 

HI  *  $kHr4p»i  *  [rJHki 

from  which  Equation  IH-2  may  be  subtracted  to  yield  the  following 
expression  for  the  errors 


in-3 


Hi  ■  [ru]H"‘l)  -  f’j 


III-4 


The  constant  stress  and  constant  strain  methods  are  now  distinguished 

(k-l) 

by  the  manner  in  which  the  are  specified.  The  Ramberg-Osgood 

stress-strain  relation,  which  we  may  write  in  the  form 

1  Vi  t  l^sr  -j^T| 

will  be  used  in  examining  both  methods. 
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Thfe  mnlinear  strains,  in  the  constant  stress  method,  are  given  by 


c(k-l) 


o  I  n(k-l)  I  n  ft(k-l) 


£3  p_  !£_ 

7EJ  Kj*  I  |^IT| 


which,  on  using  Equation  III-l,  may  be  written  in  the  expanded  form 


,<»-!) .  ^4-!  Nkl)  ° ♦  nK"'1 

J  7Ej  |Jk_1)|  la°J  oP^ 


(k-l)jn-l  _(k-l) 


cn  i  Ai 


+ - 


A  similar  expansion  can  be  constructed  for  the  e£  '  on  introducing 


Hk)]  ■  M  *Hk,i 


which  gives 


t(k)  „  _a 

J  7  E 


-oo  n.4(k-i)  |  „ 


U(k-l)|n-l  ^k) 


J  h$kJ|  [!"<* 


CT0JAJ|  aoJAj+"'"' 


Substituting  the  foregoing  expression  into  Equation  HI-4,  and  assuming 
small  errors  in  the  sense  that  and  small  steps  in  the  loading  such 

that  Aq^«q^,  so  that  terms  containing  the  products  of  these  quantities 

may  be  dropped,  we  obtain 

\A)Ur  ]f_  3jl^  °_1  U-D  M) 

^  J  L‘jJvj  |ipT|  Vj  lJ  J 


which  may  be  re-written 


|4k,|  ■  [ru]  1  -  { Sj1'1’! 

— j  j  — 


ra-5 


where  give  the  slope  of  the  stress  versus  inelastic-strain  curve  at 
the  individual  element  stress  levels. 
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Instability  develops  vhen  the  errors  in  the  step  increase  over 

those  in  the  k-1  step.  Clearly,  the  method  vould  be  expected  to  be 

(k)  (k-l) 

unstable,  if  any  of  the  inequalities  >  5^  were  satisfied 

directly.  A  more  critical  check,  however,  is  the  consider  the  entire 
set  of  as  a  vector  (in  n-dimensional  space  where  n  is  equal  to  the 

number  of  elements)  and  apply  the  condition  that  the  Eucliudean  length 
of  this  vector  does  not  increase.  Mathematically,  this  means  that  for 
stability 


where  the  critical  condition  is  defined  by  using  the  equality  sign. 

For  simplicity,  we  consider  the  case  of  infinitesimally  small  steps 
in  the  loading,  so  that  the  Aq^  terms  may  be  neglected,  and  further,  de¬ 
note 

[ru][i  '  Sj‘'1>  AJ  ’  [BU_1>] 


The  critical  condition  for  instability  then  becomes 


It  Is  now  observed  that  the  eigenvalues  of 

MT  [ 

will  all  be  positive,  hence  the  condition  that  none  of  these  eigenvalues 
be  greater  than  one  (which  becomes  the  stability  requirement) ,  can  be 
replaced  by  the  more  severe  condition  that  the  sum  of  the  eigenvalues 
not  be  greater  than  one.  This  latter  condition  can  be  assured  by  re¬ 
quiring  that  the  sum  of  the  squares  of  all  elements  of  not  be 

greater  than  one,  and  in  addition,  that  the  absolute  sum  of  any  row  or 


column  in  B,  not  be  greater  than  one.  The  critical  stress  is  then 
found  from  the  largest  of  these  sums. 


7 

In  the  case  of  the  example  truss  problem,  where  E.  =  E  =  10  psi, 

5  1 

n  =  10,  Ai  =  1.0  sq.  in.,  and  oq ^  =  oq  =  10  psi,  the  sum  of  the  squares 
approach,  viz., 


yields  =  83,800  psi.,  while  the  rows  and  columns  approach  gives 
ocr  =  82,900  psi. 

Recall  that  the  foregoing  development  has  determined  the  minimum 
value  at  which  instability  might  occur.  It  is  of  interest  to  compute 
a  critical  value  of  stress  at  which  instability  is  strongly  assured  to 
occur.  This  can  be  done  by  returning  to  the  initial  notion  that  in¬ 
stability  will  occur  if  §^)  >  .  This  amounts  to  restricting 

attention  to  the  diagonal  elements  in  the  matrix  'The  corresponding 

critical  stress  (lowest  value)  will  be  given  by  the  largest  (in  absolute 
value )  element  on  the  diagonal  in  the  B .  matrix .  This  corresponds  to 

*7  lj 

the  -  0.4l4  x  10  terra  in  the  F^  matrix,  so  that  the  critical  stress 
is  given  by 

n-1 

x(0.Wk  „  10T)  .  ! 


which  gives  ccr  =  93,800  psi.  The  lowest  stress  at  which  instability 

would  develop,  in  the  case  of  the  example  truss  problem,  would  therefore 
be  expected  to  occur  between  82,900  psi  and  93,800  psi. 

Note  that  if  only  one  structural  element  were  inelastic,  then  the 
diagonal  term  in  the  matrix  corresponding  to  this  element  would 

give  the  correct  critical  stress  by  the  latter  procedure.  Both  of  the 
foregoing  values  of  critical  stress  have  been  indicated  in  Figure  4, 
where  they  are  seen  to  correlate  with  the  experimental  (computer)  re¬ 
sults.  The  simple  approach  of  considering  only  the  diagonal  elements 
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appears  to  'be  advantageous  In  the  present  problem.  With  this  approach, 
it  is  easy  to  see  that  finite  values  of  Aq^,  which  make  the  elements  in 

the  second  diagonal  matrix  in  Equation  III- 5  less  than  one,  must  raise 
the  stress  for  instability,  vhich  explains  the  progression  of  critical 
stresses  (with  increasing  load  step  size)  appearing  in  Figure  4. 

Finally,  it  is  noted  that  certain  of  the  foregoing  results,  viz., 
the  value  of  a  based  on  one  diagonal  elements  alone,  can  be  obtained 
cr  ( v}  / v ^ 

simply  by  introducing  qi  '  +  f :  '  in  place  of  qj  '  directly  into 

^  1  (k)  ^  (k) 

Equation  HI-3,  and  regarding  the  y^  as  perturbations  on  the  q^  ' . 

The  additional  results,  such  as  the  demonstration  of  the  effects  of 
finite  Aqi,  however,  are  not  obtained  by  this  procedure. 

In  the  case  of  the  constant  strain  method,  in  addition  to  the 
error  quantities  in  the  member  loads 


we  also  introduce  error  quantities  for  the  relaxed  loads 


XV 

The  equation  defining  the  relaxed  loads,  written  for  the  generic  i 
member,  is 


>>  ,„(*)  (k)  -1 

‘Ui  3q*1f  q^ 


which  may  be  written  in  the  following  form 


If  ve  now  apply  the  condition  of  no  load  reversals,  then  the  terms 

q(k)y  jq^  |etc.,  will  all  produce  the  same  sign  (for  a  given  element) 

and  hence  may  be  cancelled  out.  Applying  the  condition  that  all  error 
terms  and  load  increments  Aq^  are  much  smaller  than  the  load  magnitude 

j  <Ljj ,  leads  to  the  result 

.(k)  .(k)  Ei  /  .(k-l)  ,(k)  A_(k)\ 

h  m  V  -  f  (k-ij  [  V  5*i  ] 


where 


|*U-i>  "‘1 

-~/v  ,  \  3n_  1  ‘V 

i/^n  *  Vi 


The  corresponding  form  of  Equation  III-4  may  now  be  written  by 
introducing  the  exact  load- reduction- increments  Aq^,  where 


which  leads  to 


,  x  qo  *(k"l)  +  p(k-l)  n  (k-l) 

( »(k))  r„  l  3  o.i  Vi  Vi  Vi 

l!i  1  N  rt.  v  1  afcn 


7  EJ  "oj  Aj  |ijw| 


Applying,  once  again,  the  smallness  requirement  on  the  and  Aq^,  yields 


TO 


where  the  absolute  value  signs  have  been  emitted  for  simplicity. 

The  check  for  the  occurrence  of  instability  may  now  be  carried  out 
in  the  same  manner  as  for  the  constant  stress  method.  Thus,  considering 
only  the  case  of  infinitesimal  load  increments,  the  corresponding  form 
of  the  matrix,  as  defined  by  Equation  III- 6,  is  found  to  be 

M  -ffeNlN 

A  numerical  check  for  the  special  case  of  the  truss  problem  shows  that 
the  critical  condition  of  the  eigenvalues  summing  to  unity  calls  for 
physically  inadmissible  values  of  A  simpler  demonstration  of 

this  property  is  provided  by  the  "direct  approach"  (i.e.,  setting 
In  this  case,  the  critical  stress  is  given  by  the 

lowest  value  corresponding  to  the  "n"  equations  obtained  by  equating 
the  diagonal  elements  of 


In  the  example  truss  problem,  where  E^  =  E,  oq^  =  oQ,  and  -  1.0, 
and  where  the  three  diagonal  elements  in  the  r>3  matrix  can  be  denoted 
by  -CjE,  where  in  turn  CK^-d,  the  foregoing  matrix  equation  reduces  to 
the  following  simple  algebraic  equation 

-y  B_  .  t!  . 

5(cr)i  +  E 

The  indicated  critical  values  are  easily  seen  to  be  -  C^E  811(1  -  (2  -  Cj)B 

both  of  which  are  physically  inadmissible  for  the  Ramberg-Osgood  stress- 
strain  relation.  Thus,  the  constant  strain  method  is  indicated  to  be 
free  of  instability  in  this  case.  The  problem  of  accuracy,  of  course, 
is  another  matter,  due  to  the  necessity  of  working  with  finite  (and  pre¬ 
ferably  large)  load  steps.  These,  apparently,  are  responsible  for  the 
slow  divergence  of  the  computer  results  shown  in  Figure  5* 
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APPENDIX  XV 

ROTATION  OP  AXES  OP  ANISOTROPY 

Let  the  x  and  y  axes  he  rotated  fran  the  axes  of  anisotropy  in  a 
positive  sense  so  that  fran  the  strain  transformation  equations  ve  get: 

d«x  =  +  m2<i«22  +  2^mdei2 

dV  -  B2d«n  +  /2 d«22  -  ^md«12 

dc^  =>  -/ md«n  t/md*^  +  2(/2  -  m2)de12 

when  /,  and  a  are  the  usual  direction,  cosines  of  the  x-axis  with  respect 
to  the  orthogonal  axes. 

Similarly  fran  the  stress  transformation  equations: 

°11  “  +  “S  ‘  ^°xy 

°22  m  m2ax+/°y  +  ^CTxy 

a,-  -  /mo  -/mo  +  (Z2  -  m2)a 
12  x  *  y  xy 

Substituting  Equation  IV-2  into  the  appropriate  expressions  in  Equation 
12  gives  a  set  of  equations  as  follows: 
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Finally  substituting  Equation  IV-3  into  Equation  IV- 1,  simplifying  and 
defining  three  new  coefficients  6,  6^  and  6g  such  that  6  = 

6^  =  (or^  +  S)  and  6g  *  (cr^  +  6)  the  appropriate  incremental  flow 
equations  become: 

dex  =  ^  ^  6i  +  m^fi2  '  4)ax  ’  ^“12  "  -  /2m26 2)ay 

-  2jt*U2^  -  n%)°xy] 

d«y  =  f5[c/46g  +  ®\  -  6)<Jy  -  (“X2  “i2m2filK 

-  2An(m261  - 

-  f[(3“lA  *  +  '•/»%)<»„ 

-  m262)crx  -  ^(m2^ 

Insertion  of  Equation  IV-2  into  11  gives  the  expression  for  effec¬ 
tive  stress: 

a2  =  [(/\  +  m462  -  6)a2  -  2^  -  ~  f^^Vy 

+  (f462  +  mkt>±  -  6)a2  +  (3o^  +  +  4/2m2a2)a2 

-  -  m262)axaxy  -  W®2^  “  /^Vxy] 

Now  suppose  we  make  an  uniaxial  stress-strain  test  for  the 
x-direction.  The  valid  expressions  becane: 


Fran  the  three  Equations  17-4,  we  obtain  equations  similar  to  Equation 
15 .  This  equation  is  augmented  under  restrictions  of  equal  plastic 
work  thus : 

a  /£■]_  =  b  C  /S£  =  +  ^*t2  -  ?  17-6 
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Assigning  the  value  of  unity  to  one  "a"  as  before  ve  obtain  the  relative 
values  of  the  others.  In  particular  ve  nov  obtain  a  value  for  (/^6, 

li  .  L 

+  o  -  6)  which  is  a  function  of  0^3  known  6114  <244 

unknown.  Therefore  a ^  may  be  determined. 


\ 
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APPENDIX  V 


INELASTIC  MATRIX  COMPUTER  PROGRAM 


A.  Program  Description 

This  is  a  brief  description  of  the  Grumman  inelastic  matrix  pro¬ 
gram  for  carrying  out  elastic-plastic-creep  analysis  -  deck  No.  45128  - 
Elastic  Unloading  (follows  Hooke's  Law  when  unloading),  including 
anisotropy . 

This  program  for  elastic-plastic  analysis  has  proved  to  be  quite 
adaptable  for  analytic  investigations.  It  has  been  modified  to  include 
an  option  for  anisotropic  analysis.  Previously,  the  program  was  modi¬ 
fied  to  use  time-hardening  theory,  deformation  theory,  inelastic  un¬ 
loading,  and  constant  stress  theory  (none  of  which  have  been  retained 
in  the  final  program).  It  was  not  necessary  to  do  major  program  re¬ 
visions  to  accommodate  these  variations.  The  program  housekeeping  is 
arranged  so  that  modifications  to  the  manipulation  of  data  will  not 
affect  the  housekeeping.  Thus  the  program  makes  a  convenient  frame¬ 
work  to  explore  various  calculation  procedures. 

The  program  is  written  in  Fortran  II  to  run  on  the  Grumman  IBM 
7094s.  These  are  2  channel  (A  and  B)  machines  with  6  drives  per 
channel ;  32,768  words  of  core  storage;  on-line  card  reader;  on-line 
printer;  and  printer  clock.  The  program  is  set  to  run  under  Fortran 
Monitor  control  which  uses  a  $J0B  card  for  identification.  Input  is 
on  logical  tape  7  (A-2),  print  output  on  logical  tape  6  (A-2)  and 
punched  output  on  logical  tape  5  (B-4)  (not  used  by  this  program). 
Logical  tape  8  (B-l)  is  used  for  storage  of  binary  output  which  is 
converted  to  BCD  print  output  in  link  6  (at  the  end  of  Job).  The  pro¬ 
gram  will  accept  an  input  data  tape  on  logical  tape  9  (A-5)  and  will 
write  a  binary  save  tape  for  restart  on  logical  tape  11  (A-6).  These 
2  auxiliary  tapes  are  optional  for  each  run  (see  description  of  the 
control  cards).  The  Grumman  IOU  subroutine,  as  well  as  the  subroutines 
for  rewinding  and  unloading  a  tape  (RUN)  and  for  moving  to  the  start 
of  a  designated  file  (FILTAP)  are  included,  as  required  in  the  program, 
in  column  binary  form. 

The  program  tape  furnished  to  Wright  Patterson  Air  Force  Base 
contains  all  the  information  needed  to  duplicate  our  analysis.  It  is 
in  the  following  sequence: 

File  1  -  a  1-card  BCD  label  tested  by  the  program 'to  distinguish 
BCD  data  tape  frcm  binary  save  tape. 

File  2  -  BCD  card  images  for  matrices  SIM  AND  SU.  See  pages  79 
and  80  for  a  description  of  the  matrices  and  their 
format . 

\ 
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Pile  3  -  BCD  card  Images  of  Fortran  program  (6  links),  each  link, 
including  binary  subroutines  previously  mentioned. 

Multiple  (  6  )  end  of  file  marks . 

We  recommend  that  all  3  files  of  this  tape  be  copied  for  use,  then 
that  file  3  be  punched  onto  Fortran  cards.  The  punching  program  must 
be  able  to  handle  mixed-mode  cards  to  accanmodate  the  short  binary  sub¬ 
routines.  This  deck  of  cards,  with  the  proper  *I.D.  or  $JOB  card,  will 
then  be  in  proper  form  for  a  Fortran  compile  and  execute,  using  the 
copied  tape  on  drive  A- 5  (logical  9)*  Tbs  program  uses  only  the  data 
fran  files  1  and  2;  it  will  not  move  into  file  3*  This  2-file  format 
for  the  data  tape  is  generated  by  the  GISMO  matrix  system  in  use  at 
Grumman  and  elsewhere. 

Each  link  of  the  program  contains  the  non-IBM  subroutines  needed 
for  operation.  Standard  input-output  subroutines  etc.  will  be  taken 
fran  the  library  tape.  As  a  point  of  information,  the  program  contains 
six  links  numbered  consecutively  fran  1  to  6. 

1)  Link  1  reads  the  first  control  card,  and  reads  all  other 
decimal  input  supplied. 

2)  Link  2  is  used  only  on  a  restart  Job.  It  reads  the  modified 
step  table,  if  provided,  and  part  of  the  binary  input  tape. 

3)  Link  3  is  used  only  on  a  restart  Job,  and  reads  the  balance 
of  the  binary  input  tape. 

k)  Link  4  is  the  processing  link.  It  does  all. the  calculation, 
print  control,  and  writing  of  binary  output  on  tape  B-l  to  be 
converted  to  BCD  print  output  by  Link  6. 

5)  Link  5  writer  a  binary  tape  for  restart,  then  transfers  to 
Link  6.  If  no  binary  tape  is  to  be  written,  exit  is  from 
Link  6. 

6)  Link  6  reads  the  binary  output  stored  on  tape  B-l  and  converts 
it  to  BCD  output  on  tape  A-3  for  printing.  When  tape  B-l  has 
been  completely  processed,  a  message  to  the  operator  indicates 
that  it  need  not  be  saved.  If,  due  to  machine  error  or 
operator  intervention,  tape  B-l  is  not  processed  into  prints 
on  A-3,  but  is  saved,  then  Link  6  can  be  used  as  a  separate 
program  using  B-l  as  a  data  tape  and  will  process  B-l  into 
prints  on  A-3.  Link  6  does  not  use  any  data  fran  COMMON.  All 
necessary  clues  are  stored  on  B-l, 

Sequencing  and  detalle  of  the  data  cards  follow.  The  symbols  used 
in  the  program  for  various  items  of  input  data  are  listed  on  page  77 
and  are  shown  on  the  sample  key-punching  sheet  page  78.  , 
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The  data  cards  are  used  in  the  following  sequence,  immediately 
after  the  *DATAcard  required  by  the  Monitor: 

1.  General  clue  card  (FORMAT  l)  containing  KLU4,  NINCLD,  HA,  KLUT, 
KLUISO,  ALPHA,  EETA,  GAMMA,  GHU,  and  a  title  or  caption. 

2.  Table  of  load  or  time  steps  desired.  Up  to  ten  cards  defining 
ten  steps  may  be  used.  Load  steps  and  time  steps  may  be  used 
in  any  sequence.  The  maximum  level  for  each  load  or  time  step 
may  be  above  or  below  the  previous  maximum  level  of  load  or 
time.  The  program  verifies  'the  algebraic  sign  of  the  increment, 
and  corrects  it  if  necessary.  Each  card  contains  four  variables 
TEMPI,  TEMP2,  TEMP3,  TEMP4  in  FORMAT  2. 

3«  Data  matrices.  These  may  be  provided  in  any  sequence.  Each 
matrix  has  a  header  card  in  FORMAT  3>  one  or  more  data  cards 
in  FORMAT  4,  and  a  blank  card  to  end  it.  The  last  input 
matrix  on  the  Monitor  input  tape  must  be  followed  by  one 
added  blank  card  (two  total)  to  trip  the  program  into  operation. 

In  the  event  that  a  job  is  running  overtime,  and  it  is  desired  to 
stop  the  program  in  a  restartable  form,  mount  a  blank  tape  on  logical 
tape  11  (drive  A-6)  and  set  Sense  Switch  6  on.  This  .will  write  the 
contents  of  memory  on  A-6  in  proper  fc—i  to  continue'  the  run  later. 

The  program  distinguishes  between  the  saved  binary  tape  and  a  decimal 
input  tape  at  starting  time.  Either  is  mounted  on  logical  tape  9 
(Drive  A-5). 

Built-in  pauses  in  the  program  are  as  follows:  ^ 

f 

1.  Pause  1.1111  to  mount  the  data  input  tape  at  the  start  of  the 
program,  if  all  matrices  are  not  on  the  Monitor  input  tape 
A-2 .  For  this  KLU4,  in  the  first  data  card,  should  be  a  "1”  to 
"4"  indicating  the  count  of  decimal  matrices  on  A-5-  If  A-5 
is  a  binary  saved  tape  from  a  preceding  run,  any  digit  (except 
zero  )acceptable  for  KLU4. 

2.  Pause  1  to  mount  a  blank  tape  on  A-6  to  receive  memory.  This 
is  reached  either  with  Sense  Switch  6,  or  with  a  card  in  the 
table  of  steps  punched  TAPE  in  columns  7-10. 

B,  Symbols  and  Format  of  the  Data  Cards 

1.  General  Clue  Card  -  FORMAT  1 

Cols .  Field  Symbol 

1  H  KLU4 


This  gives  the  number  of  input  matrices 
on  the  auxiliary  input  tape  (A-5).  If 
all  matrices  are  an  the  monitor  tape, 
leave  this  blank.  If  using  binary  re¬ 
start  tape,  use  a  digit.  Maximum  number 
of  decimal  input  matrices  on  the  auxiliary 
input  tape  is  4. 


** 
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Cols.  Field  Symbol 


2-4 

13 

NINCLD 

This  gives  the  number  of  non- increment  re¬ 
cycles  at  each  load  or  time  level.  The 
total  cycles  at  each  level  is  NINCLD  +  1. 
When  printing,  the  first,  tenth,  twentieth, 
etc . ,  and  last  cycle  of  each  printable  level 
will,  be  printed. 

5-7 

13 

NA 

This  sets  the  frame  size  for  the  problem 
to  be  handled.  NA  is  three  times  the 
number  of  nodes.  Maximum  value  is  165 
(55  nodes). 

8 

n 

KLU7 

0,  or  blank,  prints  5  preselected  matrices 
on  cycles  indicated  by  the  step  table; 

1  prints  ell  matrices  on  cycles  indicated 
by  the  step  table.  2  prints  5  preselected 
matrices  on  all  cycles;  3-9  print  all 
matrices  on  all  cycles. 

9 

11 

KLUISO 

0  indicates  an  isotropic  run 

1  indicates  an  anisotropic  run 

10 

IX 

Not  used 

n-20 

E10.3 

ALPHA 

A  variable  defined  by  the  creep-strain 
equation. 

21-30 

E10.3 

BETA 

A  variable  defined  by  the  creep-strain 
equation. 

31-35 

F5.2 

GAMMA 

A  variable  defined  by  the  creep-strain 
equation. 

36-40 

F5.2 

GNU 

Poisson’s  ratio  "nu" 

41-50 

10X 

Not  used 

51-80 

5A6 

TA-TE 

Any  30  characters  of  alpha-numeric  text 
to  be  printed  as  a  heading  for  identifies- 

tion  purposes. 

2.  Table  of  Steps  (limited  to  ten  entries)  -  FORMAT  2 


Cols • 
1-6 


Not  used 


LOAD  indicates  a  load  step,  TIME  indicates  a 
time  step,  TAPE  indicates  write  memory  on  a 
save-tape  on  drive  A- 6  (logical  #11 )  then  exit. 


7-10 


Field 

6X 

A4 
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Cola.  Field 

FIHS  or  blank  Indicates  end  of  table  (this  may 
be  the  lltfc  card  in  the  table). 

H-20  E10.6  Upper  limit  of  step  in  pounds  or  hours 

21-30  E10.1  Interval  or  increment  for  calculation 

31-40  ELO.l  Interval  or  increment  for  print  output.  Prints 

are  generated  on  tape  A- 3  for  the  current  cycle 
,  when  the  current  load  or  time  level  is  an  inte¬ 

gral  multiple  of  the  print  interval.  If  the 
print  interval  is  left  blank,  no  prints  are 
generated  for  cycles  in  this  step.  If  print 
interval  is  very  small  compared  to  the'  current 
level,  then  numeric  problems  sometimes  occur 
in  the  print  control  subroutine  (OUTPUT),  and 
it  may  be  necessary  to  re-run  with  the  every- 
cycle  print  control  "2"  or  "3”  for  KLU7  punched 
in  the  first  control  card. 

41-80  Ignored 

3.  Data  Matrices  -  Header  Card  -  FORMAT  3 


Cols . 

Field 

1-6 

6x 

Not  used 

7-io 

4x 

Not  used.  We  use  the  letters  MTRX  for  com¬ 
patibility  with  the  GISMO  Matrix  System,  which 
reads  and  writes  matrices  in  this  format. 

•11 

VC 

Not  used 

12-17 

A  6 

This  is  the  identification  name  for  the  input 

matrices  and  must  correspond  exactly  with  one 
of.  the  following  names: 


bbbSIM  Matrix  of  stresses  for  applied  loads 
maximum  size  165x1 

bbbSLT  Matrix  of  stresses  for  member  strains 
maximum  size  165x165 

bTKIGN  Table  of  stress  values  llxl 

bfTSPSN  Table  of  strain  values  llxl 

These  two  matrices  define  the  stress- 
strain  curve  as  a  series  of  chords. 
The  data  is  entered  in  this  format 
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merely  to  conform  to  the  format  of  the 
83M-SU  matrices  vhich  were  generated 
using  the  GISMO  Matrix  System.  Note 
that  the  first  value  in  both  TSIGN  an 
TEPSN  must  be  zero  to  avoid  upsetting 
the  interpolation  procedure. 

TALF  12  o^2  ) 

TALF  a  ) 

J  23  <  Basic  anisotropic  parameters 

<  un 

TALF  31 

TALF  44  j 


18 

IX 

Not  used 

19-21 

13 

Number  of  rows  in  this  matrix 

22-24 

13 

Number  of  columns  in  this  matrix 

25-32 

8X 

Not  used 

33 

11 

The  digit  3 

34-80 

Not  used 

Data  Matrices  - 

Data  Cards  -  FORMAT  4 

Cols. 

Field 

1 

IX 

Not  used 

2-4 

13 

Row  index  for  the  first  element 

5-7 

13 

Column  index  for  the  first  element 

8-23 

EL6.8 

The  first  element  on  this  card 

24 

IX 

Not  used 

25-27 

13 

Row  index  for  the  second  element 

28-30 

13 

Column  index  for  the  second  element 

31-46 

E16.8 

The  second  element  on  this  card 

47 

IX 

Not  used 

48-50 

13 

Row  index  for  the  third  element 

51-53 

13 

Column  index  far  the  third  element 
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Field 


* 

Cols. 

% 

54-69-  i’,'  E16.8  The  third  element  on  this  card 

V  *’ 

70-80  Hot  used 

The  last  card  of  a  matrix  must  be  completely  blank  (tested  in  Col. 
2-4).  i  The  last  matrix  on  the  Monitor  input  tape  7  (drive  A-2)  must  be 
followed  by  one  added  blank  card  (two  total)  to  trip  the  program  into 
operation. 

The  input  matrices  on  the  Monitor  tape  may  be  in  any  sequence  as 
long  as  each  matrix  starts  with  a  header  card,  has  all  its  data  cards 
next,  and  ends  with  a  blank  card. 

C.  Anisotropic  Parameter  Matrices 

If  a  run  is  indicated  as  isotropic  in  the  first  control  card 
(KLUISO  in  Col.  9  is  zero  or  blank)  the  program  will  read  in  the 
anisotropic  parameters,  if  provided;  then  it  will  replace  them  with  the 
built-in  parameters  for  the  isotropic  condition.  If  a  run  is  indicated 
as  anisotropic  in  the  first  control  card  (KLUISO  Is  1  to  9)  >  then  one 
each  of  the  anisotropic  paraneter  matrices  must  be  provided  or  the  pro¬ 
gram  will  terminate  on  an  error. 

D.  Restart  Procedure 

If  the  run  being  set  up  is  a  restart,  the  input  deck  can  be  in 
several  forms.  The  first  control  card  must  have  a  digit  in  KLU4  so 
that  the  program  will  read  tape  A-5;  HINCLD  and  KLU7  are  read  from  this 
card.  The  other  factors  are  carried  from  the  previous  run.  The  con¬ 
tinuation  of  an  Isotropic  run  will  always  be  isotropic,  and  conversely, 
regardless  of  the  clue  provided  (KLUISO). 

The  table  of  steps  may  be  read  in  again  (modified)  if  the  previous 
run  is  stopped  with  sense  switch  6,  or  it  may  be  retained  and  continued 
from  the  previous  run.  However,  if  the  previous  run  is  stored  on  tape 
A- 6  by  using  a  TAPE  card  in  the  step  table,  then  a  new  step  table  must 
be  read  in.  "™" 

In  either  case,  the  last  data  card  on  a  restart  job  must  be  FINS 
or  blank  in  columns  7-10.  This  means  a  restart  data  deck  will  have  a 
minimum  of  two  cards  (clue  card  and  a  blank),  or  a  maximum  of  12  cards 
(clue  card,  10  step  cards  and  a  FINS  or  blank). 

E.  Time  Estimates 


For  time  estimates,  allow  3.5  minutes  to  compile  the  Fortran, 

1.5  minutes  to  read  the  input  tape,  2.0  seconds  per  cycle  printed,  and 
100  to  110  cycles  of  calculation  per  minute.  Each  printed  cycle  writes 
approximately  3  feet  of  print  tape  (55  node  problems). 
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MATRIX  ANALYSIS  FOR  ISOTRoPlC  UK  ANiSOTjROPIC  INELASTIC  STRUCTURES 
IN  THE  PLASTIC  AND  CREEP  RtGlMt 


THIS  PROGRAM  WRITTEN  FOR  a  IK  FuKLfc  CONTRACT  AF  33(6151-2260 
*  *  *  TABlE  UF  SYMBOLS  USED  •  *  * 

*  *  ARRAYS  *  * 

ALL212  -  MATRIX  OF  AN1SUTKUPIC  PARAMETERS  FOR  EACH  NODE 

AL1223  -  MATRIX  OF  ANISOTKuPiw  PARAMETERS  FOR  EACH  NUDE 

AL12T1  -  MATRIX  OF  ANISOTRUPli  PARAMETERS  FOR  EACH  NUOE 

ALFA44  -  MATRIX  OF  ANlSUTRuPlC.  PARAMETERS_FOR  EACH  NUDE 

OELEPK  -  MATRIX  OF  NUDE' STRAIN  CHANGES' I  DELTA  EPS  IllONl 

FOR  THE  X  AND  Y  DIRECTIONS,  AND  THE  SHEAR  STRAIN  CHANGE 
DELEPN  -  MATRIX  OF  EFFECTIVE  INELASTIC  STRAIN  CHANGES 
EPBARN  -  MATRIX  OF  EFFECTIVE  iNELAstIC  STRAINS  (EPSILON  3AR  SUB  NT 
FOR  K-TH  CYCLE 

EPBARP  -  MATRIX  OF  EFFECTIVE  INELASTIC  STRAINS  (tPSILON  BAR  SUB  N> 
FOR  K-l  CYCLE  ( THE  PRECEDING  CYCLE)' 

EPCNK  -  MATRIX  OF  EFFECTIVE  uKEcP  STRAINS  FOR  THIS  CYCLE 

EPCNP  -  MATRIX  OF  EFFECIIVE  CREEP  STRAINS  FOR  PREVIOUS  CYCLE 

EPSUK  -  MATRIX  OF  NUDE  PLASIIL  STRAINS  (EPSILON  SUB  U)  FUR  THE 
X,  Y  AND  SHEAR  STRAINS,  FUR  THE  K-TH  CYCLE 
KPMTM  -  TABLE  OF  CONTROL  CLUES  J LOAD  STEP,  TJME  STEP*  WRITE  SAVE 
TAPE)  READ— I N  CONTROL  INFORMATION'  USED  WITH  PMTM 

PMTM  -  TABLE  OF  LOAO  INCREMENTS,  PRINT  CONTROL  INCREMENTS  ANO 
UPPER  LOAD  LEVEL  PER  SIEP  (CONTROL  INFORMATION) 

SGBARM  -  MATRIX  OF  EFFECTIVE  NuDE  STRESSES  FOR  LAST  CYCLE  THAT 
SHOWED  AN  INCREASE  Al  A  PARTICULAR  NODE 
SGBARN  -  MATRIX  OF  EFFECTIVE  NuDE  STRESSES  FOR  CURJtfcNT  CYCLE 

SGBARP  -  MATRIX  OF  EFFECTIVE  NuDE  STRESSES  FUR  PREVIOUS  CYCLE 

SIGUK  -  MATRIX  OF  NODE  STRESSES  (SIGMA  SUB  U)  FOR  K-TH  CYCLE 
SIJ  -  MATRIX  OF  STRESSES  FOR  MEMBER  STRAINS 

SIM  -  MATRIX  OF  STRESSES  FUR  APPLIED  LOADS 

TALFI2  -  MATRIX  UF  ANISOTROPIC  PARAMETERS  (INPUT  DATA} 

TALF23  -  MATRIX.  OF  ANISOTROPIC  PARAMETERS  tINPUT  DATA) 

TALF3I  -  MATRIX  OF  ANISOIROPIC  PARAMETERS  (INPUT  DATA) 

TACF44  -  MATRIX  OF  ANISOIROPIC  PARAMETERS  (INPUT  DATA) 

TEFSTN  -  MATRIX  OF  TOfAL  EFFECTIVE  STRAINS 

TEPSN  -  TABLE  OF  STRAIN  VALUES  DEFINING  THE  STRESS-STRAIN  CURVE 
TIMK1  -  MATRIX  OF  REFERENCE  CREEP  TIMES  FOR  ALL  NOOES 

TOTEPS  -  MATRJX  OF  TOJAL  NUDE  STRAINS _ _ 

TSIGN  -  TABLE  OF  STRESS  LEVtLS  DEFINING  THE  STRESS-STRAP'  CURVE 

•  *  VARIABLES  AND  CLUES  *  * 

ALPHA  -  PARAMETER  USED  IN  THE  CRfcEP  EQUATIONS 

BETA  -  PARAMETER  USED  IN  THE  CREEP  EQUATIONS 

E  -  MODULUS  OF  ELASTICITY _ 

GAMMA  -  PARAMETER  USED  IN  THfc  CREEP  EQUATIONS 

GNU  -  POISSONS  RATIO 

K  -  THE  CYCLE  COUNTER 

KERRSW  -  CLUE  USED  FOR  TEMPORARY  INDICATOR  BETWEEN  LINKS 

KLU4  -  CLUE  INDICATING  TOTAL  COUNT  OF  MATRICES  ON  AUXILIARY  TAPE 

KLUS  ~  CLUE  INDICATING  MATRICES  STILL  NOT  REAP  FROM  AUXIL.  TAPE 


GOO  l 
UU02 
0003 


0005 
COOS 
0007 
0008 
0009 
0010 
0011 
CD  12 
0013 
0014 
0015 
001b 
0017 
00  IB 
0019 
0020 
0021 
092? 
0023 
0024 
0925 
0026 
0027 
002  8 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0937 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
00S6 
0057 
0058 
0059 
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KLU6  -  CLUE  FOR  PRINT  LUNfRUL  I  WHICH  CYCLES)  0060 

KLU8  -  CLUE  FOR  PRINT  CUtfTnUL  (WHICH  MATRICES)  0061 

KLUISO  -  CLUE  FUR  ISOTRuPlc  uK  ANISOTROPIC  RUN  (INPUT)  00b2 

KSET  -  CLUE  FOR  THE  LOAJ  UK  Tint  LEVEL  CURRENTLY  IN  USE  0063 

NA  -  3  TIMES  THE  NOOt  lUUNT  (INPUT)  0069 

NC  -  NUMBER  OF  NOOES  f)0b5 

NINCLU  -  NUMBER  OF  NUN-  iNwKEMtNT  CYCLES  AT  EACH  LOAD  LEVEL  (INPUT)  00  ' 

PM  -CURRENT  LOAD  LEVEL  .•  00o7 

SHRMOO  -  SHEAR  MODULUS  OObB 

TIME  -  CURRENT  TIME  LEVEL  0069 

0070 

OUTPUT  -  SUBROUTINE  TO  CunThul"  PRINTING  AT  VARYING  LOAD  LEVELS  0071 

0072 

'  0073 

ELASTIC  UNLOAOING  ( FOLLOWS  .HUuKcS  LAM  WHEN  UNLOADING)  0079 

COMMON  TEFSTN,  TOTEPS  •  007V 

COMMON  KLU9,  KLU5,  KLiJo ,  KLU8  >  NA,  NC,  K  0U76 

COMMON  KERRSM,  NINCLO,  KSET,  PM,  E,  GNU,  SHRMUU  007 f 

COMMON  ALPHA,  BETA,  GAMMA,  TIME,  KLUISO  0078 

COMMON  PMTM,  KPMTM,  >IM,  SIJ,  TSIGN,  TEPSN  0079 

COMMON  TALFI2,  TALF23,  IALF31,  TALF99  0080 

COMMON  AL1212,  AL1223,  ALI231,  ALFA99  0081 

01MENSI0N  TALF12I II ) »  T  ALP^s( LI) , I ALF31 ( II) , TALF99( 11)  0082 

DIMENSION  AL 1 2121 55  I ,AL 12231 55) , ALL 231 1  55) ,ALF 999(55)  0083 

EQUIVALENCE  ( T IMKl, TEFSTn) ,( Tul gPS,0£LEPK )  __  0089 

THESE  FOUR  ARRAYS  ARE  EQUIVALENCES)  JO  SAVE  CQR~E  SPACE.  0085 

01  MENS  ICN  TIMKK55)  , TEFS1 N( 55 1 , TuIEPSC  165  )  ,DELEPK(  165 )  0086 

0 1  MENS  ION  PMTM<L0,3),KPMTMU0),SIM( 165),SIJ(l65,165),TSIGN(ll)  0087 

DIMENSION  TEPSN(il)  0088 

DIMENSION  BL(5I,I0(9),NR(3), NU( 31 , EL ( 3 )  0089 

DIMENSION  PSTEP <  3 )  0090 

EQUIVALENCE  ( ID(l), ADI), (10(2) , A02) , ( ID ( 3 ) , A03) , ( 10(9 ) , A09 ) , ( 10(5)  0091 

1,AU5I ,( ID(6),A06) ,( (0(7)  , AU7 ) , ( lu( 8) • AD8 ) , ( 10(9) , AD9)  0092 

ADI  *  606060623199  0093 

AD2  «  60606062 J 191  0099 

A03  -  606362312795  0095 

A09  *  606325976295  _  _  _  0096 

A05  ■  632193260102  0097 

A06  »  632193260203  0098 

AD7  *  632193260301  0099 

AOS  *  632193260909  0100 

AD9  *  606060606060  0101 

BL(l)  *  939621296060  _  _  _ _  0102 

BL ( 2)  »  633199256660  .  0103 

BL( 3)  ■  632197256060  0109 

BL( 9)  ■  263195626060  0105 

BL ( 5)  ■  606060606060  0106 

KERRSM  *  l  0107 

INTAPE  *__7  _  _  .  0108 

KLU5  •  6  0109 

KLU6  »  0  0110 

KLU8  *  10000  0111 

KLUALF  ■  0  0112 

PSTEP(l)  -  6.0  0113 

PSTEP(2)  »  0.0  _  _  ' _  0119 

PSTEP (3)  «  6.0  *  0115 
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00  51  J22»JUI0  _  _  __  0116 

51  KPMTMI J22)  »  0  0117 

LROW  *  0  0118 

READ  INPUT  TAPE  7, 1 ,KLU4,N1NClU ,NA,KLU7, KLUISO, ALPHA, 8E T A, GAMMA,  0119 

1GNU,  TA«TB*TC»TO»T£  0120 

WRITE  OUTPUT  TAPE  6, 2 1 , TA, Td, IC , T U, TE  0121 

KLU7  1  INPUT)  -  0  TO  PRINT  PKt-SEl.EC  TED  MATRICES  UN  SELECTED  CYCLES  0122 

KLU7  1  INPUT)  *  l  TO  PRINT  ALL  MATRICES  ON  SELECTED  CYCLES  0123 

KLU7  (INPUT)  *  2  TO  PRINT  PKt-SELECTEO  MATRICES  UN  ALL  CYCLES  0126 

KLU7  (INPUT)  -  3-9  TO  PRINT  ALL  MATRICES  ON  ALL  CYCLES  0125 

I F  ( KLU7 1384,384, 381  0126 

381  IF(KLU7-1) 383,383,382  0127 

382  KLU6  *  1  __  ; _  _  0128 

KLU6  (OUTPUT)  *  l  TO  PRINT  t  ACH  CY'cLE,  0  TO  PRINT  SELECTED  CYCLES  0129 

IFIKLU7-2) 384,384,383  0130 

383  KLU8  «  0  0131 

KLU8  (OUTPUT)  *  0  TO  PR1NI  ALL  MATRICES  0132 

KLU8  (OUTPUT)  *  10000  TU  PRINT  itLECTED  MATRICES  0133 

384  CONTINUE  0134 

1 F ( KLU4 ) 302 ,302,301  0135 

301  PRINT  17  0136 

PAUSE  11111  0137 

TEST  INPUT  TAPE  TO  DISTINGUISH  GISMU  BCD  FROM  INELASTIC  PROG.  SAVE  0138 
REAO  INPUT  TAPE  9,31,LTAPt  0139 

REWIND  9  _  _  0140 

SUBROUTINE  FILTAP  POSITIUNi  A  TAPt  AT  THE  FIRST  RECORD  OF  ANY  FILE  0141 
CALL  F ILTAP ( 9,2 )  0142 

IF(LT APE— 21)343, 343 ,344  0143 

344  CONTINUE  0144 

SAVED  TAPE  FROM  PREVIOUS  RUN  UF  INELASTIC  PLATF  0145 

THIS  RUN  IS  A  CONTINUATION  0146 

CALL  CHAIN  (2,3)  0147 

343  CONTINUE  0148 

GISMO  FORMAT  BCO  TAPE  0149 

THIS  RUN  IS  A  NEW  ONE  0150 

KLU5  *  KLU4  0151 

INTAPE  *  9  _ _  _  0152 

302  CONTINUE  0153 

KLUISO  *  0  FOR  ISOTROPIC  RUN  0154 

KLUISO  -  1  FOR  ANISOTROPIC  RUN  0155 

WRITE  OUTPUT  TAPE  6,9  0156 

IF( KLUISO) 380,380,379  0157 

379  CONTINUE  __  0158 

WRITE  OUTPUT  TAPE  6,10  0159 

380  CONTINUE  0160 

52  READ  INPUT  TAPE  7,2 , TEMPI, ltMP2 ,TEMP3,TEMP4  0161 

LROW  =*  LROW  ♦  l  0162 

00  53  J22«l,5  0163 

B  IF(TEMPl*(-BL(  J221_)  )53,54153 _ __  _  0164 

53  CONTINUE  0165 

C  BAO  CONTROL  CARO  0166 

WRITE  OUTPUT  TAPE  6, 12, LROW  0167 

WRITE  OUTPUT  TAPE  6,2.TEMP1,TEMP2,IEMP3, TEMP4  0168 

CALL  EXIT  0169 

54  GO  TO (55, 55 j 55,60,60)  ,322 _ 0170 

55  IF(LROW-IO) 56,56, S3  0171 
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KPMTMt  N  I  *  l  FOR  *  LOAk  >UP 

KPMTMt  N  )  *  2  POR  A  I  lMt  t  a o c 

KPHTHC  N  )  *  3  TO  DUMP  HU«uKY  l«iu  ^  iAVt  TAP€ 

PMTM  <N,1)  «  UPPER  LIMIT  uf  **1l^  CALCULATION 

z:  ir,: :  \~  hkss"  k  sas1*-™ 

KPMTMILROMI  *  J22 

PMTM< LRGW, 1 J*  TEMP2 

tEM(>5  *  TEHP2-PSTLPI J22> 

PMTHILRCM* 21*  SIONFI TEMP At lENPal 
PSTEPl J22)  »  TEMP2 


TEHP5  *  TEHP2-PS I tr l 

PMTM(LR0M,2I*  SIONFI TEMPAt ItHPsi 
PSTEPt J22 J  »  TEMP2 
IF (PMTM I  LRON»2I  157*63*57 
GO  TO<995,995,57>  ,  J22 
CONTINUE 

PHTMILR0N«3I*  TEMP4 

GO  TO  52 

CONTINUE 

DO  361  J23  *  1.10 

IMKPMTMU23)  1361,361,362 
1FIKPMTHI J231-41 363,36i,iol 
J24  *  KPMTMI J23I 

GO  TO  (364,366,368) , J24  Mthi n 

WRITE  OUTPUT  TAPE  6,22 .PHIMI J24 .21 # PMTMU23U 1 

IF(PMTM(J23. 31)361, 361, 3o5 

WRITE  OUTPUT  TAPE  6 ,23  ,PMIH(  J<43* 4*  — 

WRITE  OUTPUT  TAPE  6 ,24  , PHIMI J23 *2» tPHTHI J23,  1 > 

IFIPNTNIJ23, 311361, 361, 3o7  ,  • 

'  WRITE  OUTPUT  TAPE  6,26,PMIrt( 323,3) 

GO  TO  361 

l  WRITE  UUTPUT  TAPE  6,27  . - 

L  CONTINUE  „„  4  , 

WRITE  OUTPUT  TAPE  6,28,NlNwL0 
WRITE  OUTPUT  TAPE  6,25,  ALPHA, Bel  A.uAHMA 
IF l GNU 1 86, 37 ,87 
i  GNU  *  .3 

?  CONTINUE  .  .  . 

IF  I NA 1999,91, 92 

1  NA  *  165 

2  IFINA-165193, 93,999 

3  CONTINUE 

6  IF(NA-3*INA/3) 1999,97,999 

7  NC  *  NA/3  .  -  -  ~ . " 

IF ( KLU4 1 304 , 304, 303 

3  IF  I KLU4— 4) 306,306,996 
6  CONTINUE 
INTAPE  *  9 

>8  READ  INPUT5TAPE  INTAPE,  3, NAMc, NRUwJ, , NCOL S , NF ORH 

00  110  121*1, »  ,, ,  ,  . 

IF! NAME- ID! 121 1 1 110,31 l.AlO 

L°  8A0T INPUT  -  MATRIX  NAME  No!  AuCtPTAdLE 

11  GO  TO  (321, 322,323,324,323,320^327,^28, till j 121 

21  00  331  ii  •lflN* 


0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 

0181 

0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0190 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 

0211 

0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 
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331  SIMt Ill  »  0.0 _ 

WRITE  OUTPUT  TAPE  6,26,N*Me,NRUw:»,N(.aLS,  INTAPE 

GU  TO  111 

322  00  332  11  -1.165 
00  332  12  >1.165 

332  S1JII1. 12)  «  0.0 

GO  TO  111 _  _  __  _ 

323  00*333  II  »  1,11 

333  TSIGN(U)  -  0.0 
GO  TO  ill 

324  DO  334  11  >  1,11 

334  TEPSN(Il)  >  0.0 

GO  TO  111  _ _ 

325  DO  335  II  *  l.ll 

335  TALF12( II )>  0.0 
KLUALF  -  KLUALF  4-1 
GO  TO  111 

326  DO  336  II  >  1,11 

336  TALF231 1 1 ) »  0.0 _ _ 

KLUALF  »  KLUALF  ♦  2 

GO  TO  111 

327  DO  337  11  >  1,11 

337  TALF31I 1 1 ) >  O.C 
KLUALF  *  KLUALF  ♦  4 

GO  TO  IU  _ _ _ 

328  00  338  11  >  1,11 

338  TALF44I Ill-  0.0 
KLUALF  »  KLUALF  +8 
GO  TO  111 

111  GO  TO  (112,112, 112, 112, 112, 112,112,1 12, 1501, 121 

112  READ  INPUT  TAPE  INTAPE ,4, (NR 1 122 ) ,N0( 122 ) ,EL(  122 ), 122-1  .3 ) 
IF ( Nftl 1 ) >996,109,113 

109  I F ( KLU5 1308,308, 306 
308  INTAPE  >  7 
GO  TO  108 
304  GO  TO  308 

113  GO  TO  (121. 122, 125, 126, 201, 202, 203. 204, 150), 121 _ 

C  READ  IN  ARRAY  SIGMA-IN 

121  MROW  »  NR(1)  . 

WRITE  OUTPUT  f APE  6,4,(NRU22),N0(122),EL( I22),I22>1,3) 

SIM  ( MROW )  *  ELI  1 )  _  _ 

IF (EL (2 >1127, 128, 127 

127  MROW  >  NR(  2 ) _ _ _ : _ _ _ 

SIM  (MROW)  >  EL ( 2 ) 

128  IF (ELI 3) >129.112. 129  _  _ 

129  MROW  >  NR (3) 

SIM  (MROW)  «  EU3I  _  .  ... 

GO  TO  112 

C  READ  IN  ARRAY  SIGMA- 1 J  _ _ ; _ 

122  MROW  >  NR( 1 ) 

MCOL  -  N0(  1)  _  _ _  _  .  . .  . 

SIJ  (MROW, MCOL)  >  EL(  l) 

IF  (EL  (2)  1130, 131^130  _ _  .  ..... 

130  MROW  ■  MR(2I 

MCOL  •  N0(2) _ _ 

SIJ  (MROW, MCOL)  »  EL(2) 


0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 

0261 

0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 

0278 

0279 

0280 

0281 

0282 

0283 


131  IF l EL t  3) I  132,  112,  132  0284 

132  MROW  *  NR(3)  0285 

MCQL  *  N0(3»  0286 

S1J  (MROW.MCQLl  =  EL  (31  0267 

GO  TO  112  0288 

C  READ  IN  ARRAY  TS1GN  <  TABlc  Jt-  ^ lv>RA  BAR  Nl  0289 

125  MROW  *  NR  ID  0290 

TS1GNIMR0M )  *  ELI1I  0291 

1F(EL(21 >139, 140,139  0292 

139  MROW  »  NR  t  2  )  0293 

TSIGN(MROW 1  -  EL  t  2 1  0294 

140  IF(EL(3) 1 141. 112.141  0295 

141  MROW  NR(  3)  0296 

TS1GMMRUW1  *  EL131  0297 

GO  TO  112  0798 

C  READ  IN  ARRAY  TEPSN  (TArtLt  Ul-  tPSlLUN  BAR  N)  0299 

126  MROW  >  Nftlll  0300 

TEPSNtMRUWl  «  ELI  1 )  0301 

IF ( EL ( 2 ) 1 142. 143. 142  0302 

142  MROW  -  NR(2)  0303 

TEPSNCMROWI  *  EL(2)  0304 

143  IF(EL(3)I 144. 112. 144  0305 

144  MROW  >  NR ( 3  )  0306 

TEPSNCMROWI  ■  EH3I  0307 

GO  TO  112  0308 

C  READ  IN  ALPHA  TABLES  0309 

201  MROW  =*  NRCll  0310 

TALF121MRUW1  *  ELlll  0311 

1F(EL(2) 1221.222.221  0312 

221  MROW  *  NR( 2  )  0313 

TALF121MR0W)  *  ELI21  0314 

222  1F(EL(31)223,112,223  0315 

223  MROW  *  NR( 3)  0316 

TALF121MR0W1  *  EL( 3  I  0317 

GO  TO  112  0318 

202  MROW  *  NR(11  0319 

TALF23CMROWJ  *  EL(1»  _ _ _  _  '  0320 

IF(EL(2) 1224,225.224  0321 

224  MROW  *  NR ( 2  1  0322 

TALF23CMR0W)  ■  EL(21  ,  0323 

225  IF(EL(3) 1226.112,226  0324 

226  MROW  -  NR( 3)  0325 

TALF231 MROW 1  *  EL(31  . . . . .  _  _  0326 

GO  TO  112  0327 

203  MROW  *  NR (11  0328 

TALF311MR0W1  *  ELUI  0329 

IFIELI2) 1227.228,227  0330 

227  MROW  *  NR( 2)  0331 

TALF3KMR0WI  *  EL(2)  0332 

228  IFIEL131 1229,112,229  0333 

229  MROW  *  NR(31  0334 

TALF3KMR0W1  *  EL(31  0335 

GO  TO  112  0336 

204  MROW  ■  NRUl  0337 

TALF44IMR0W1  -  EH11  _ _ _  _  _  0338 

1F(EL(21 1230,231,230  *  0339 
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230  MROW  »  NRJ21  _ _ _  _ _  0340 

TALF44!  MROW  1  *"  EL  1 2 1  '  "*  ~ 0341 

231  IF ( EL ( 31 )232«  1 12  t  232  0342 

232  MROW  *  NR ( 3 )  0343 

TALF44! MROW )  *  ELI3)  0344 

GO  TO  112  0345 

150  CONTINUE  _  _ _ _ _  0346 

1 F ( KLU I  SOI  1 65  * 165,161  “  ***'  ”  ~  ~  0347 

161  IF{KLUALF-15) 16?. 167, 162  0348 

162  WRITE  OUTPUT  TAPE  6,32,KLUALF  ’  0349 

IF!KLU4)164,164, 163  0350 

C  SUBROUTINE  RUN  REWINDS  AND  UNLuAUS  THE  DESIGNATED  TAPE  0351 

163  CALL  RUN! 91  0352 

164  CALL  EXIT  ~  ~  *'*  . *  0353 

165  00  166  I  *  1.11  0354 

TALF12I 11*  0.5  0355 

TALF23I 1 1  *  0.5  0356 

TALF31I Il»  0.5  0357 

166  TALF44I  1 1  *  1.0  _  _ _ _  0358 

167  CONTINUE’  '  . .  ■  *  “  0359 

E  *  TSIGNI21/TEPSNI2I  0360 

81  SHRMOO  *  E/I2.6I l.+GNUI 1  0361 

WRITE  OUTPUT  TAPE  6,5,E,SHKMU0,GNU  0362 

WRITE  OUTPUT  TAPE  6,6  0363 

DO  149  Il*lxll  _  _ _ _  __  _  0364 

149  WRITE  OUTPUT  TAPE  6.7,11, f SluNi 111,1 EPSN III 1, TALF12I lil,TALF23(ll)  0365 

l.TALF31(Ill,T  ALF44I  1 1 1  0366 

K  -  0  0367 

PM  *  0.0  0368 

IFIKLU41152.152.209  0369 

C  SUBROUTINE  RUN  REWINOS  AND  UNLUADS_THE  DESIGNATED  TAPE  0370 

209  CALL  RUN(  9 1  '  "  *  . . .  0371 

151  PRINT  19  0372 

152  CALL  CHAIN ( 4,31  0373 

994  WRITE  OUTPUT  TAPE  6,18,KEU<,  0374 

CALL  EXIT  0375 

995  WRITE  OUTPUT  TAPE  6,16,LkOW _ _ _  _  0376 

CALL  Exit  ~  0377 

996  CONTINUE  0378 

WRITE  OUTPUT  TAPE  6,14,NAMt  *  0379 

WRITE  OUTPUT  TAPE  6,4,  I NRU221  ,NOi  122 1  ,EL  1 122 1 ,  I22=>1 , 3  1  0380 

CALL  EXIT  0381 

99  7  CONTINUE  _  _  _ , _  _  0382 

WRI  E  OUTPUT  TAPE  6,13V  NAMt  -  •  ‘  0383 

CALi  EXIT  0384 

999  WRITE  OUTPUT  TAPE  6,ll,NA  0385 

CALL  EXIT  0386 

1  FORMAT! 11,213,211, IX, 2E10. 3, 2F5. 2, 10X.5A6I  0387 

2  FORMAT  1 6X ,  A4.  E10. 6,  E  IQ.l.ElO.ll _ _  0388 

3  FORMAT !  I IX,  A6,  IX 72  1 3,8X~»  1 1 J  " .  0389 

4  FORMAT! 3( IX, 213, E 16. 81  1  0390 

5  FORMAT ( 26H1  MODULUS  OF  ELASi ICltv  ■  ,F11.0,4H  PSI ,6X, 16HSHEAR  MODU  0391 

1LUS  »  , FI  1. 0 ,4H  PSI .6X.5HNU  »  ,F6.3J  0392 

6  FORMAT I41H0  TABLE  OF  VALUES  FUK  STRESS-STRAIN  CURVE  //  0393 

15X.29HP0INT  STRESS  LEVEL  SlKAlN,7X,8HALPHA  12,7X,8HALPHA  23, 7X  0394 
1,8HALPHA  31, 7X.8HALPHA  44/16X ,3HPSl ,9X , 7HIN./IN.// 1  0395 
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7  PUKMAr<oXtl3»2X,Fll.2o!3A,Fl2.dl)  0396 

9  FORMAT  (  //  79H  THIS  t  SUTruP  l  t  run  USES  ELASTIC  UNLUAUlNG  (HOOKES  L  0397 
1AM)  mITH  STRAIN  HAROCN l No)  Q398 

10  FORMAT) 1HF,5X,2HAN>  0399 

11  FORMAT ( 10H  ERROR  NA=I4)  0400 

12  FORMAT { 2bH  ERROR-  INCREMcn!  uAru  nu.,13,5H  N.G. )  0401 

11  FORMAT  (  14H  ERROR-MATRIX  ,Ao)  (H02 

14  FORMAT ( 17H  EAROR-NFG. INDEX  A»i  0401 

16  FORMAT  I 33H  ERROR-NO  IN  Ter Vau-1  nLhEMENT  CARO, I  3)  0404 

17  FORMAT ( S IHO  PAUSE  lllll  lu  MuOM  INPUT  DATA  lAPt  UN  DRIVE  A-S//// I  0405 

lb  FORMAT ( 7H  ERROR- 15, 18H  InPuT  MaTrICeS  NG)  040o 

19  FORMAT  1 U9H  OEMUUNT  ANO  SAVt  TaPe  A-3.  OAT  A  HAS  SEEN  READ  INTO  MEMO  0407 

1RY.  THE  PROGRAM  CONTINUE©  IU  RUN.//)  0408 

20  FQRMAT(//5X»7HMATR1X  , Ao, 1A, 1 i, ort  KuMS  X  ,I3,l4H  COLUMNS  FROM  TAPt  0409 

I  .12)  0410 

21  FORMAT! 1H1,29X,5A6)  0411 

22  FORMAT (5X, 16HL0AD  INCReMenIS  ,FV.2,llH  POUNDS  TO  ,F10.2,7H  PJUNOS )  0412 

23  FURMATI 1H+ , 6 IX , l 9HPR 1 N T  uOIPUl  EVeKY  ,F9.2,7H  POUNDS)  0413 

24  FURMATI 5X, 16HT IME  INCREMENTS  ,r4.4,llH  HOURS  TO  ,Fl0.4,hH  HOURS)  0414 

25  FORMAT!  /  9H  ALPHA  =*  , E 10. J, 5X , 7HuE I  A  *  , t 10 .3 , 3X , bHGAMMu  *  ,E10.  3)  0415 

'26  FORMAT! 1H*,61X,19HPRINT  OUTPUT  EVERY  ,F9.4,oH  HOURS)  0416 

27  FUAMAT(5X,35HSTORE  MEMORY  UN  IAPE  A-6,  THEN  EXIT)  0417 

28  FURMAT(//5X,I3,48H  NUN- InURcMenT  LYULES  AT  EACH  LOAD  OR  TIME  LEVEL  0418 

II  0419 

31  FORMAT  UOX,  12)  0420 

32  F0RMAT15X, 1 2 , 65H  ALPHA  Ta8LeS  WERE  READ  IN  UF  4  REQUIRED,  CHECK  YU  0421 

1UR  INPUT  CAROS. I  0422 

END(l,l,0,0,0,0,l,l,0,l,0,U,0,u,0) 

*  CHAIN! 2,3)  ,  0424 

*  LIST8  0425 
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o  o  o  o  eo  to  as  so  a> 


<151282  GRUMMAN  AIRCRAFT  ENoIrtEfcHlNG  CuRP.  DECK  NU.  45128  LINK  2 


C4S1232  MATRIX  ANALYSIS  OF  INELASTIC  PcATE  -  LINK  2  -  WITH  CREEP  0427 

C  THIS  LINK  READS  IN  A  SAVtUBiNARY  TAPE  -  FIRST  PART  0428 

.  COMMON  TEFSTN,  TOTEPS  0429 

COMMON  KLU4 ,  KLU5 ,  KlJ6,  KcU8,  NA,  NC ,  K  0430 

COMMON  KERRSW,  NINCLD,  KStl,  PM.  E,  GNU,  SHRMOJ  0431 

COMMON  ALPHA,  SETA,  uAMMA ,  TIME,  KLOISO  0*32 

COMMON  PMTM,  KPMTM,  S1M»_  S1J,  TSIGN,  TEPSN  0433 

COMMON  TALF12,  TALF23,  TAL.F31 ,  IALF44  ‘  0434 

COMMON  AL1212,  AL1223,  AL1231,  A*.FA44  0435 

DIMENSION  TALF12(ll),TALF23(ll),IAt.r31(m,TALF44tll)  043b 

DIMENSION  AL  1.212  (  55  ) ,  At  1223 1  55  ),AL  1231(55),  ALF  A44  ( 55  )  0437 

COMMON  SIGUK,  EPSUK  0438 

EQUIVALENCE  I TIMK1 , TEF STNl , I  TUI tPS.OELEPK)  0439 

C  THESE  FOUR  ARRAIfS  ARE  EQUlvALtNCtu'lO  SAVE  CORE  SPACE.  0440 

DIMENSION  TIMK1I55) , TEPSINI 55) , IufEPSI 165) ,UELEPK 1 165 )  0441 

DIMENSION  PMTMI 10,3) , KPMTM l lU)«SiM(165),SIJ(l66,l65),TSIGN(ll)  0442 

DIMENSION  TEPSNI 11) ,S IGUK1 lu5 ) ,EPSUK( 165)  0443 

DIMENSION  8LI5)  0444 

DIMENSION  PSTEPI 3)  0445 

READ  TAPE  9,KLU4,KLU5,  ”'lA,  L  8 » N A , NC  *  K  *  K  E  RR  S W  j  ~LC»KS£T,PM,fc,'  0446 
IGNU, SHRMOO, ALPHA, 3ETA, GAMMA, I1ME,KLUIS0  0447 

KLU6  *  1  WILL  PRINT  EVERY  cYCle,  U  WILL  PRINT  ONLY  SELECTED  CYCLES  0443 
KLU8  *  0  WILL  PRINT  ALL  MATRICES,  NUN-ZERO  PRINTS  SELECTED  MATRICES  0449 
NINCLO  (NON-INCREMENT  CYCLES  AT  EACH  LOAD  LEVEL)  ACCEPTED  FROM  0450 

INPUT  CONTROL  CARO  _ _ 0451 

WRITE  OUTPUT  TAPE  6^9  ~  '  “  .  0452 

IFIKLUt SO 1380, 38 0,379  0453 

379  CONTINUE  0454 

WRITE  OUTPUT  TAPE  6,10  0455 

380  CUNTINUE  0456 

NCLU  a  1  _  _  __  _  _  _  _  0457 

KERRSW  «  2  ' .  ‘  ”  ~ . .  045  5 

BL (13  *  434621246060  0459 

BLI2)  »  633144256060  0460 

8L( 3)  »  632147256060  0461 

8L( 4)  »  263145626060  0462 

8U5J  *  606060606060 _ 0463 

00  51  142  *  1,10  ~  0464 

51  KPMTM! 142)  *  0  0465 

PSTEPI l )  *  0.0  0466 

P$TEP(2)  »  0.0  0467 

PSTEPI 3)  =  0.0  0468 

LROW  *  0  . _  0469 

52  READ  INPUT  TAPE  7.«2, TEMPI,  IEMP2,  TfcMP3,TEMP4  . . .  ~  0470 

LRUW  a  LRUW  +  1  0471 

DO  53  J22al,5  ‘  “  04? 2 

8  IF(TEMPl*(-BL(J22)) 153,54,63  0473 

53  CONTINUE  0474 

C  BAD  CONTROL  CARO _ 0475 

WRlfE  OUTPUT  TAPE  6,12,LKUw"  '  . .  0476 

WRITE  OUTPUT  TAPE  6,2,TEMPl,TfcMP2,TEMP3,TEMP4  0477 

CALL  EXIT  0478 

54  GO  TO I  55,55,55,60,601 ,  J22  0479 

55  I F (  LROW— 10)56,56,53  . .  0480 

C  KPMTM!  N  )  a  1  FOR  A  LOAD  STEP _ 0481 
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O  J  J  J  u 


KPMTHI  N  »  «  2  FUR  A  TIMc  ileP  0*82 

KPMTM(  N  )  *  3  TO  OUHP  MtMoKY  tNlu  A  SAVE  TAPE  0*83 

PHTH  (N* l »  «  UPPfcR  LIMIT  uF  SltP  046* 

PHTH  IN«2)  *  INTERVAL  (  IiYuKcMcNT I  rUK  CALCULATION  0*85 

PMTM  ( N • 3 1  *  INTERVAL  IlHuKtHfcNl)  FuR  PRINT  OUTPUT  0*86 

56  KPMTMILKQW)  *  J22  0*6/ 

NCLU  a  2  0*68 

KSET  *  1  0*69 

PMTMl  LRUW, 1 )*  TEMP2  0h90 

TEMP5  *  TE  MP2-PS IEPI J22 )  0*91 

PHTMI LROWt 2 )*  SIGNFI TEMPJ, IEMP5I  0*92 

PSTEPI J22 1  *  TEMP2  0*93 

IF(PMTM(LRUW»2) )57»63*57  0*9* 

63  GO  TO  (995,995,57),J22  0*95 

57  CUNTINUE  0*96 

PHTMI LROWt  3 ) ■  TEMP*  0*97 

GO  TO  52  0*98 

60  CONTINUE  0*99 

GO  TO  161,62) ,NCLU  0600 

62  CONTINUE  *  0501 

READ  TAPE  9,JUNK  0502 

READ  TAPE  9, JUNK  0503 

GO  TO  6*  050* 

61  CONTINUE  0505 

READ  TAPE  9,PMTM  0506 

READ  TAPE  9.KPMTM  0507 

6*  CUNTINUE  0508 

WRITE  OUTPUT  TAPE  6,1  0509 

DO  361  J23  -  1,10  0510 

IFIKPMTMIJ23) 1361,361,362  0511 

3o2  IF ( KPMTHI J23 )-*) 363,361,601  0512 

363  J2*  »  KPMTHI J23 )  0513 

GO  TO  (36*, 366, 368), J2*  051* 

36*  WRITE  OUTPUT  TAPE  6,22tPM!M(J2J,2)»PMTM(J23,l)  0515 

1F1PMTHIJ23, 31)361,361 ,3o5  0516 

365  WRITE  OUTPUT  TAPE  6 , 23 , PHTMI J2J ,3)  0517 

GO  TO  361  0518 

366  WRITE  OUTPUT  TAPE  6 ,2* , PHTMI J23,2) ,PMTMf J23, 1 )  0519 

IFIPMTMIJ23, 3)1361, 361, 3o7  0620 

367  WRITE  OUTPUT  TAPE  6, 25, PHTMI J2i,3)  0521 

GO  TO  361  0522 

368  WRITE  OUTPUT  TAPE  6,27  0523 

361  CONTINUE  052* 

WRITE  OUTPUT  TAPE  6, 26,N1NCLU  ’  "  . . .  0525 

READ  TAPE  9, SIM  *  0526 

READ  TAPE  9, SI J  0527 

READ  TAPE  9.TSIGN  0528 

READ  TAPE  9.TEPSN  0529 

REAO  TAPE  9,TALF12  0530 

READ  TAPE  9,TALF23  “  .  .  0531 

REAO  TAPE  9  ,TALF3 1  0532 

REAO  TAPE  9,TALF**  0533 

REAO  TAPE  9«AL12] 2  053* 

REAO  TAPE  9.AL1223  0535 

REAO  TAPE  9,AL123J. _ 0536 

READ  f  APE  9,  ALFA**  *  ~  '  "  0537' 


READ  TAPE  9.SIGUK  _ _  _  0538 

RE AO  TAPE  9 ,EP SUK  '  ""  "  0539 

CALL  CHAIN  (3,3)  0540 

995  WRITE  OUTPUT  TAPE  6,16, LKUW  0541 

CALL  EXIT  0542 

1  FORMAT ( IH0,29X,3 7HC0NT 1NUAI I  UN  kU«  -  INELASTIC  ANALYSIS)  0543 

2  F0RMAT(6X,A4,E10.6,E10.1 , tl O.i) _  _  0544 

9  FORMAT ( // 79  H  THIS  ISOTKUPIC  RUN  USES  ELASTIC  UNLOADING  (HOOKES  L  0545 

1 AW)  WITH  STRAIN  HARDENING  05*6 

10  FORMAT! 1H*,5X,2HAN)  0547 

12  FORMAT ( 26H  ERROR-  INCREMENT  CARO  NU.,I3,5H  N.G.)  0548 

16  FORMAT ( 33H  ERROR-NU  l NTEHVAL- INCREMENT  CARO, I  3)  0549 

22  FORMAT ( 5X* 16HL0A0  INCREMENTS  ,F9.2,11H  POUNDS  TO  ,F10.2,7H  POUNOSJ  0550 

23  FORMATt lH»,6lX,i9HPRINT  OuTPUl  EVERY  , F9. 2 ,7H  POUNDS)  .  0551 

24  FORMAT ( 5X, 16HTIME  INCREMENTS  ,F9.4,11H  HOURS  TO  ,F10.4,6H  HOURS)  0552 

25  FORMAT  (  1H+,61X,1 9HPR I  NT  OUTPUT  EVtkY  ,F9.4,6H  HOURS)  0553 

26  FORMAT ( //5X, I 3,48H  NON-INCREMENT  CYCLES  AT  EACH  LOAO  OR  TIME  LEVEL  0554 

l!  0555 

27  FORMAT!  5X_,35HST0RE  MEMORY  uN  TAPE  A-6,  THEN  EXIT) _ _  0556 

ENO( 1,1, 0,0 ,0,6, 1,1,6, 1,0,0, U, 0,0) 


* 

* 


CHA IN ( 3,3) 
LIST8 


0558 

0559 


451283  GRUMMAN  AIRCRAFT  ENGlNfctRlNo  cuhP 


DEC*.  NJ.  45128  LINK  3 


C451283  MATRIX  ANALYSIS  OF  INllmSIiC  PlAI £  LINK  3  -  KITH  CREEP  0561 

C  THIS  IS  THE  SECOND  HALF  uE  uLU  lINk  2  0562 

C  THIS  LINK  READS  IN  A  SAVtu  uiNAKY  InPE  -  SECOND  PART  0563 

COMMON  TEFSTN,  TOTEPS  0564 

COMMON  KLU4,  KLU5,  Kuuo,  KLU8,  NA ,  NC ,  K  0565 

COMMON  KEKRSw *  NINCLO,  KjE I ,  PM,  E,  GNU,  SHKMoJ  0566 

COMMON  ALPHA,  OETA,  uAMMa ,  TIME,  KLUISO  0>67 

COMMON  PMTM,  KPMTM,  sin,  SlJ,  TSIGN,  TEPSN  0568 

COMMON  T ALP  12  ,  TALF23,  TALEjI,  TALF44  0564 

COMMON  AL1212,  AH223,  AUdii,  ALfA-,4  0870 

DIMENSION  TALF12I  11),TALE*3( m ,IAlE3K ll),TALF44Cll>  0571 

DIMENSION  ALI212(55), ALl2z3i  33) , Ac  12 31 (  55) *  ALE  A44(  85)  0572 

COMMON  SIGUK ,  EPSUK  0573 

COMMON  SG8ARN,  SGBAKP ,  SudAKh,  EPdAKN,  EPbARP ,  OELEPN  0574 

COMMON  EPCNK,  EPCNP  0575 

EQUIVALENCE  1 TIMK1, TEFSTN),! loTtPS.OELEPK)  0576 

c  these  four  arrays  are  eqoiyalenced  to  save  cure  space.  0577 

D IMCNSI  ON  T1MK1I  551  ,  TEFSlnOs)  ,IuIcPS(  165)  ,OELEPK(  1651  057tt 

D I  MENS  I  UN  PMTMI 10,31 .KPMIMI io ) , SI M 1 165 ) , S 1 J I  165, 165 ) , TS IGN 111  I  0574 

DIMENSION  TEPSN ( III , S IGUM L65 ) ,tF SUM  165 ) , SGSARNI 55 ) , SGHARP I  85 1  0560 

DIMENSION  SOB  ARM  (551  ,  EPb«R.\|{  55  J  ,  cPbARP  (  55 )  ,  DELEON  (  55 )  ,  EPCNK  (  55)  0581 

DIMENSION  EPCNP ( 55 )  0582 

READ  TAPE  4 , SGBAKN  0583 

READ  TAPE  9.SG0ARP  0584 

READ  TAPE  9, SG8ARM  0585 

READ  TAPE  9.EP8ARN  0586 

READ  TAPE  9, EPBARP  0587 

READ  TAPE  9.DELEPN  0588 

READ  TAPE  9, EPCNK  0589 

READ  TAPE  9, EPCNP  0590 

C  SUBROUTINE  RUN  REWINDS  AuD  UNLOADS  THE  DESIGNATED  TAPE  0591 

CALL  RUN( 9 )  0592 

PRINT  31  0593 

WRITE  OUTPUT  TAPE  6, 5 ,E , SHKMUD.wNU  0544 

WRITE  UUTPUT  TAPE  6,6  0545 

DO  149  11*1,11  0596 

149  WRITE  OUTPUT  TAPE  6, 7, 11, !S1 oft till*  TEPSN! Ill, TALE  12 (11), TALF2 3(11)  0597 

1, TALF31 (Il),TALF44IIll  0548 

CALL  CHAIN  (4,3)  0599 

5  FORMAT ( 26H0  MODULUS  OF  EtAiTIGilY  *  ,F11.0,4H  PSI ,6X, 16HSHFAR  MOUU  C600 

1LUS  *  »Ell.0,4H  PSI , 6X, 5HNJ  *  ,F6.3)  0601 

6  FORMAT  1 41H0  TABLE  OF  VALDES  fUk  STRESS-STRAIN  CURVE  //  0602 

15X , 29HP01NT  STRESS  LEVEL  SIKaiN , 7X, 8HALPHA  12 , 7X, 8HALPHA  23, 7X  0603 

1 , 8H ALPHA  3  1 , 7 X,  8HAL  PHA  4<t/ioX, 3HPS1  , 9X, 7H l N. / 1 N.// )  0604 

7  FORMAT! 6X,I3»2X*F11.2,5(3X*F12.8) )  0605 

31  FORMAT ( 1 18H0I F  NECESSARY  lu  STuP  IH1S  RUN  BECAUSE  IT  EXCEEDS  THE  T  0606 

II ME  ESTIMATE,  PUT  A  RING  It*  THE  sAVc  TAPE  ON  A-5  AND  CHANGC  IT  /  0607 

256H  TO  DRIVE  A-6  TO  UPDATE  IT,  IrttN  PUT  SENSE  SWITCH  6  ON.  //)  0608 

END( 1, 1,0, 0,0, 0,1,1, 0,1, 0,0,U,D,U) 

*  0610 

*  SAVE  TAPE  a-i  UNLESS  ON-LINE  PRINT  SAYS  IT  HAS  BEEN  PROCESSED  0611 

*  0612 

*  PUT  SENSE  SWITCH  6  ON  TO  END  THE  RUN  IF  IT  EXCEEDS  THE  TIME  EST1M.  0613 

*  0614 

*  IF  PAUSE  l  OCCURS  (SENSE  SWITCH  o  UR  INTERNAL  CONTROL) ,  MOUNT  A  0615 

*  BLANK  TAPE  ON  A-6.  THIS  lAPE  MILL  HAVE  RESTART  DATA  WRITTEN  ON  IT,  0616 

*  AND  MUST  BE  SAVED.  0617 

*  CHA IN ( 4, 3 )  '  0618 

*  LISTS  0619 


451284  GRUMMAN  AIRCRAFT  CNGlNttKiNU  CdKP.  DECK  NU.  4512d  LINK  4  _ 

C451264  MATRIX  ANALYSIS  OF  INEcAal IV  PLATE  -  LINK  4  -  WITH  CREEP  0621 

C  THIS  LINK  JOES  THE  CALCULmIIUN  ai*U  WRITES  PRINT  OUTPUT  0622 

COMMON  TEFSTN,  T0T6PS  0623 

COMMON  KLU4 «  KLU5,  KlUo,  KlU8,  NA,  NC ,  K  P624 

CUMMON  KERRSW*  NINCLD,  KSfc),  PM,  E,  GNU,  SHRMUD  0626 

COMMON  ALPHA,  BETA,  uAMMA ,  TIME,  KLUISO  0o26 

CUMMON  OMTM,  KPMTM,  ilMt  SiJ,  TSIGN,  TfcPSN  0627 

COMMON  TALFI2,  TALF23,  TALI- 31,  TALF44  0628 

COMMON  ALI212,  AL1223,  Al1231,  ALFA44  0629 

DIMENSION  TALF12I II  I, T ALF2 III , T ALE 31 1  1 1 )  ,TALF44 1 1 1 1  OoiO 

DIMENSION  AL12I2<55),ALl22.»(53),rtLlZ31l55>,ALF444(55)  Oo3I 

COMMON  SIGUK,  EPSUK  0632 

COMMON  SGSAKN ,  SG8ARP ,  SodARM,  tPttARNi.  EP8ARP,  OELEPN  0633 

COMMON  EPCNK ,  EPCNP  0634 

EQUIVALENCE  ( TIMKI , TEESTfi)  , 1  TUT tPS.DELEPK )  0635 

C  THESE  FOUR  ARRAYS  ARE  EQdl VALfcNCtD  TO  SAVE  CORE  SPACE.  0636 

C  THEY  ARE  NOT  THE  SAME,  BUT  ARfc  Nul  NEEDED  AT  THE  SAME  TIME  0637 

C  AND  00  NOT  CARRY  FRUM  CYCLE  Tu  CYCLE.  0636 

0IMENS1UN  T IMK1 I  551  ,TEFS1M53),  Turj;P3I165  )  .OELEPKI 1651  0j39 

DIMENSION  PMTMI  10 , 3  I  ,  KPM  TMUU )  ,  S I M 1 1 65  )  ,  S  I J  (  165 , 165  ) ,  T  S  IGN  III )  0a<*0 

DIMENSION  TEPSNt  11)  ,  SIGUK  I  L  03)  tEPsuKl  1651 ,  SG8ARNI  55 »  ,  SGBARPI  551  0o4l 

DIMENSION  SG8ARM (551 *  EPBARNI  55) ,ePbARP(55) .DELhPNI 55) *EPCNKl35)  064 2 

DIMENSION  EPCNPI55J  0643 

KSET  *  KSET  0644 

KERRSW  *  K£RR$W  _  _  _  0645 

GO  TO  1418, 419), KERRSH  0r>46 

C  INITIALIZE  WORK  AREAS  0o47 

418  DO  102  II  -  1,165  0648 

EPSUK  III)  a  0.0  0649 

102  SIGUK III)  *  0.0  0650 

DU  103  II  *  1,55  _  ‘  0651 

OELEPNI  ID  *  0.0  0652 

SG8ARNI Ill  *  0.0  0653 

SG8ARPI 1 1 )  *0.0  0654 

SGBARMU1)  *  0.0  0655 

EPBARP 111)  *  0.0  0o56 

EPCNKI  ID  *  0.0  _  _  0657 

EPCNP I  II)  *  0.0  0658 

AL1212I  ID*  2.*TALF12(2)  0659 

AL1223I  ID*  TALF12I2)  ♦  TALF23C2)  0660 

AL123 HID*  TALF12I2)  ♦  IACF3U2)  0661 

■ALFA44I I  II*  TALF44I 2 )  0662 

103  EPBARNI  ID  *0.0  .  _  _  0663 

KSET  *  l  0664 

TIME  *  0.0  0665 

419  CONTINUE  0666 

REWIND  8  0667 

KSET  *  KSET  0668 

KERRSW  *  KERRSW  _  _ _ _  _  _  •_  0669 

NINCPM  *  NINCLD  ♦  2  0670 

NINTUT  *  NINCLD  ♦  l  0671 

151  K  »  K  ♦  1  0672 

KLU2  *  1  0673 

C  IF  KLU2  *  0,  THE  CYCLE  UPtHAIluNS  WILL  BE  PRINTED  0674 

C  KLU6  *  1  WILL  PRINT  EVERY  CYCLE.  0  WILL  PRINT  ONLY  SELECTED  CYCLES  0676 


105 


on  on 


ALL  MAUICEa.  NUN-2ERU  PHIMIS  SELECIED  UN£S 


KLU8  *  0  WILL  PRINT 

1F(KLU6>248,248,249 
249  KL  J2  3  0 
248  CUNT  I NUE 

IF(K- 1)270,270,271 
27 u  KLU2  *  0 
271  CONTINUE 

GO  TO  (416,417), KERRSW 
417  KLU2  3  0 
KERRSW  3  l 

KSET  IS  THE  ROW  OF  KPMTH  u*  PMTM  dcING  USED  (CURRENT  LOAD  OH 
TIME  LEVEL) 

300  IKKSET-ll)301,319,3l9 

301  IF(KPMTM(KSET))302.302.3U3 

302  KSl I  *  KSET  ♦  l 
GO  TO  300 

303  IF(KPMTM(KSET)-4)304, 302, 302  fv( . - 

VARIABLE  (J24)  INDICATES  A  LUAu  gYl.lE  (D*  A  TIME  CYCLE  (2), 
OR  WRITE  MEMORY  ON  A  SAVE  TAPE  (3) 

304  J24  *  KPMTM(KSET) 

KSET  3  KSET 

KERRSW  3  KERRSW 

IF(PMTM(KS£T,2) ) 342 ,342,341 

341  GU  TO  l 305* 315,423) *J24 

342  GU  TO  ( 345.355,423) ,J24 

305  IF (PHTM( KSET, 1 ) - ( PM*PMTM(KStT ,2) 1 1306,307, 308 

343  1F(PMTM(KSET,  l)-(  PM*PMTMUScT»2) )  (308,307,306  , 

307  IF(NINCLD)309,309,308 

309  KLU2  3  0 

303  OELPM  3  PMTM( KSET ,2 ) 

GU  TO  154 

306  1F(NINCPM-2)316,316,154 

316  KLU2  3  0 

GO  TO  302  u 

315  IF(PMTM(KSET,l)-(TIME*PNTM(*SeT,2)))306,317,31H 

355  IF( PMTM (KSET, i>-(TIME+PMTM(KS£f,/>) >318,317,306 

317  IF(NINCL0I314, 314,318 

314  KLU2  *  0  > 

318  OELT 1M=  PMTM( KSET ,2 ) 

GO  TO  330 

319  IF(KLU2)205,205,320 

320  KL02  3  0 

KSET  3  KSET  -  1 
IF(KSET)205,205,321 

321  IF(KPMTM(KSET))320, 320,3/2 

322  J24  =  KPMTM(KSET) 

GO  TO  ( 308,318) ,J24 

3  Vi  °NT  I  N'JE 
1.,  3TINUE 
33.  ONTINUE 

IFININCLO >501,501, 510 

510  I F  ( K- 1 )  501 »  50 1 » 5 1 1 

511  NINCPM  3  NINCPM  -  1 
IFINlNCPM-31513,515,515 

513  CONTINUE 


0676 
0677 
0678 
0679 
0680 
0681 
0682 
0683 
0684 
0685 
0o86 
06b  7 
0688 
0684 
0690 
0691 
0692 
0693 
0694 
0695 
0696 
0697 
0698 
0699 
0  700 
0701 
0702 
0703 
0704 
0705 
0  706 
0707 
0  708 
0709 
0710 
0711 
0712 
0713 
0714 
0715 
0716 
0717 
0718 
0719 
0720 
0721 
0722 
0723 
0724 
0  725 
0726 
0727 
0728 
0729 
0730 
0731 


106 


o  n  o  non 


GO  TO  (501,512) , NINCPM  0732 

NINCSW  «  l  MEANS  THIS  IS  A  Kt-uYCi-t  STEP  (NO  LOAD  OR  TIME  INCREASE)  0733 

NIN-S*  *  2  MEANS  THIS  IS  lAS!  RE-wYLLE  AT  THIS  LOAD  UR  TIME  LEVEL  0734 

NINES*  *  3  MEANS  THIS  IS  A  LUAU  OK  I  HE  INCREMENT  STEP  0735 

512  NINCSW  *  2  0736 

GO  TO  514  0737 

515  NINCSW  -  1  073R 

GO  TO  514  0739 

501  GO  TO  (502,503), J24  0740 

502  PM  *  PM  ♦  OELPM  0741 

GO  TO  504  0742 

503  TIME  =*  TIME  ♦  0ELT1M  0743 

504  NINCPM  *  N1NCL0  *2  0744 

NINCSW  *  3  0745 

514  CONTINUE  0746 

IF (NINCLO) 524, 524,525  0747 

525  NINKCY  *  NINCLO  ♦  3  -  NINCPM  0748 

GO  TO  (523,524,524) ,N1NCS*  0749 

523  NINKLU  «  NINKCY  -  (NINKCY/ lUi*lU  0750 

1F( NINKLU)328,524,328  0751 

NINKCY  IS  COUNT  OF  TOTAL  CYCLcS  AT  IHIS  LOAD  OR  TIME  LEVEL  0752 

NINKLU  *  0  WILL  PRINT  EAcH  TEnTh  CYCLE,  IF  THIS  LOAD  OR  0753 

TIME  LEVEL  IS  TO  8E  PRINTEU  0754 

524  GO  TO  ( 324,325) ,J24  0755 

324  CALL  0UTPUT(PM,PMTM(KSET,3)tKLyU_  0756 

GO  TO  326  0757 

325  CALL  OUTPUT (TlME,PMTM(KSt(,3) , KLU1)  0758 

326  IF( KLU2 >327,328,327  0759 

327  KLU2  =  KLUl  0760 

328  CONTINUE  0761 

IF(TIME)333,333, 332  .  _  0762 

332  GO  TU (3 33,3371,324  0763 

C  COMPUTE  REFERENCE  CREEP  T irtt  FuR  ALL  NODES  0764 

337  00  331  13  -  l,NC  0765 

EPCNPI 1 31  •*  EPCNKII3)  0766 

FACNUM  »  ABSF( EPCNPI 13))  0767 

EXPON  »  BETA*ABSF(SGBARN(i3l)  0768 

FAC0N1  »(  2.71828183**EXP0N)-1.  0769 

TIMK1 (13)  ■  (FACNUM/I ALPHA*FACDNl) >**(1./GAMMAI  0770 

GO  TO  (336, 336,334), NINCs*  0771 

336  CONTINUE  0772 

EPCNKII3)  *  ( ALPHA* (TIMK 1U3)  ) **GAMMA ) *FACDN1  0773 

GO  TO  331  0774 

334  CONTINUE  0775 

EPCNKI 13)  *  ( ALPHA* { TIMKK  131+OtLT 1M)**GAMMA)*FAC0N1  0776 

331  CONTINUE  0777 

333  CONTINUE  0778 

8  156  FA  «  214747433125  0779 

KLU9  »  l  0780 

1F(KLU2>250, 250,251  0781 

250  CONTINUE  0782 

C  LABEL  -  APPLIED  LOAD  0783 

B  FB  *  246043462124  0784 

B  FC  *  606060606060  0785 

KLUSIZ  ■  2  _  _ _ _ _  0786 

WRITE  TAPE  B, KLU9, KLUS IZ, F A, FB, PC, PH, TIME  0787 


107 


no  «(Bon  o  o  ccoaoon  o  o  ct  w  o>  a  w  o  o 


251  CONTINUE  0706 

KLU9  «  2  0769 

IF(KLU2*KLU8)252,252,253  0750 

252  CONTINUE  0791 

LABEL  *  REF.  CREEP  TIME  079 2 

FA  *  512526336023  0793 

FB  *  5125254760o3  0794 

FC  =  314425606060  0795 

WRITE  TAPE  8.KLU9,  Nt,rM,FB,FL,U  IMK1  ( J1 I , Ji=l,NC»  0796 

FA  *  252626336023  079B 

FB  *  512525476062  0799 

FC  *  635121314560  0800 

WRITE  TAPE  8.KLU9,  NC ,Fm ,Fb ,FL t i EPCNK  I J 1 ) , J 1= 1 ,NC )  0801 

253  CONTINUE  0902 

SAVE  PRECEDING  CYCLE  VALUes  ul-  fcPoARN  0803 

DO  152  II  »  1,55  0804 

152  EPuARP (ID  *  EPBARNI  1  1 )  0805 

CALCULATE  NUDE  STRESSES  -  MAIR1A  sli»UK  -  FRAME  SUE  165  X  1  0806 

DU  661  15  »  1 tNA  0807 

S1GUKI  1 51  =»  S  IM(  151  *PM  0808 

DU  861  16  =  l.NA  0809 

SIGUMI5I  *  SIGUKII5)  ♦  SUI 15,  lo)*tPSUKI  16)  0810 

861  CONTINUE  0811 

IFIKLU2J254, 254,255  0812 

254  CUNTINUE  0613 

LABEL  *  NODE  STRESSES  0814 

FA  *  454624256062  0815 

FB  *  635125626225  0816 

FC  *  626060606060  »  0817 

WRITE  TAPE  8.KLU9,  NA ,FA«F?t FC , t SIGUK  I J 1 ) . Jl*l ,NA 1  •  C818 

255  CONTINUE  0819 

CALCULATE  MAGNITUDE  ANU  Slw  UF  EFFECTIVE  STRESS  AT  EACH  NUDE  0820 

CALCULATE  EFFECTIVE  STRESS  -  MATRIX  SGBAKN  -  FRAME  SUE  55  X  l  0821 

DO  166  17  *  1 »NC  0822 

SGBARPI 17)  »  SGBARNI 17)  0823 

M3  *  3*17  0824 

M32  »  M3— 2  0825 

M3 l  *  M3-1  0826 

SGBARNI I7)>SQRTF(AL1231U/)*Sl6UR(M32)**2-AL1212II7)*SIGUK(M32)«  0827 

IS IGUK I M31 I *AL 12231 17 )* SIGUK I M 31 J**2*3.*ALFA44 ( 17) *SIGUK (M3) **2 1  0828 

166  CONTINUE  0829 

IFIKLU2)256,256,257  0830 

256  CONTINUE  0831 

LABEL  =>  EFF.  STRESSES  0832 

FA  *  252626336062  0833 

FB  »  635125626225  0834 

FC  *  626060606060  0835 

WRITE  TAPE  8.KLU9,  NC fFA,F8>FC , 1 SGBARNC J1 )  ,J  1=1,  NC  I  0836 

257  CONTINUE  0837 

CALCULATE  EFFECTIVE  INELASTIC  STRAIN  FOR  EACH  NUDE  -  INTERPOLATE  0838 

IN  TABLE  (TSIGN  VS.  TEPSN)  0839 

DO  181  18  »  1 *NC  0340 

C  SGBARP  IS  EFFECTIVE  STRESS  UF  PKEV1UUS  CYCLE  0841 

1FISGBARNI 181-SGBARPI 18)1411, 401,401  0842 

C  EFFECTIVE  STRESS  IS  ABOVE  PREVIOUS  LEVEL  0843 

C  SGBARM  IS  EFFECTIVE  STRESS  UF  LAs  I  CYCLE  TO  SHOW  AN  INCREASE  0844 
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401  IFISGBAKNI  18)-SG8ARMI  l 8 > )4i 1 ,4U2» 402  0H45 

C  EFFECTIVE  STRESS  IS  ABUVt  KNtt  UF  PREVIUUS  DRUP-UFF,  IF  ANY  9346 

402  IF  l  SGBAKP I  1  8)  -SOB ARM!  I  8  )  J  tO J  , 40-> , 403  0847 

40 J  SG8ARMII8)  *  SCBARN (18)  0848 

OU  171  19  =  1,11  0049 

ESUPRK  * (  SG8ARN (  1 8  )  /  E  )  ♦  tPttARPlib)  9860 

IFIESUPRK-TEPSNI 19) >173,1/2,1/1  Od51 

171  CONTINUE  ‘  '  085? 

GO  TO  998  0858 

172  8 ARSON  =  TSIGNII9)  0354 

IF  I SGBARNI 18I-TSIGNI2I)  1/0,177,17/  0855 

177  AL1212I 18) =  2  ,*TALF12I  19)  035o 

AL1223I  18)-=  T  AL  F 1 2 1  I9)>  TALF23I  191  9857 

AL1231I I  d )  =  TALF12I  19)*  lALFsillV)  0d58 

ALFA44II8)=  TALF44I 19)  0359 

178  CONTINUE  0360 

GO  TO  174  0861 

173  KKK2  *  19  0862 

KKK 1  *  19  -  1  0863 

STNRAT  =  I ESUPRK- TEPSNIKKKI) )/i TtPSN I KKK2)- TEPSNIKKKI) )  0364 

BARSGN  =  TS1GNIKKK1 )♦! I SiGN l KKK2 >-TS I GNI KKK1 ) )*SINRAT  0365 

IFI SGBARNI 18) -TSIGNI 21  )  17o,175,l/5  0866 

175  CONTINUE  0867 

ALFA12  =  TALF 12  I KKKl )  -M  TACP 12 1 KKK2 )-TALF 12 1 KKKl ) )*STNRAT  0868 

ALFA23  =  TALF23(KKK1  )*ITALF23URK^i-JALF23(KKRl)  )*STNRAT  0869 

ALFA31  =  TALF31 IKKKl )♦( I ALF31IKKK2)-TALF31 (KKKl )  >*STNRAT  0870 

ALFA44I I8)=TALF44IKKK1)*IIaLF441 KKK2 ) -TALF 44 l KKKl ) I *STNKAT  0871 

AL 121 2 1  18 )  =*  2  .*ALFA  12  0372 

AL1223I 1 8) 3  ALFA12  ♦  ALFA28  0873 

AL1231 1  18)=  ALFA12  ♦  ALFAJ1  0874 

176  CONTINUE  0878 

GO  TO  174  0876 

174  EPBARNII8)  =  ESUPRK  -  BARSGN  /  E  0877 

C  CALCULATE  TOTAL  EFFECT  IVc  STRAIN  -  MATRIX  TEFSTN  -  FRAME  SIZE  55X1  087b 

TEFSTNI 18 )  *  ESUPRK  0879 

415  TEFSTNI  18)  *'  TEFSTNI  18)  ♦  fcPCNKtlB)  0880 

CALCULATE  EFFECTIVE  STRAIN  CHANGES  -  MATRIX  OELEPN  -  FRAME  55  X  1  0881 

CALCULATE  INCREMENTAL  EFFECT  I VE  INELASTIC  STRAIN  0382 

OELEPNI 18 )  =  EPBARNII8)  -  EPBARPI 18)  0883 

GO  TO  1420,425), J24  0884 

425  OELEPNI 18)  =  OELEPNI 18)  +  EPCNKI 1 81  -  EPCNPI 18)  0885 

420  CONTINUE  0886 

GU  TO  181  0887 

DROP-OFF  OF  EFFECTIVE  STKtbS  0888 

OR  STILL  BELOM  THE  KNEE  UF  PREVIOUS  OROP-OFF  0889 

411  EPBARNI 18)  =  EPBARPI 18)  0390 

TEFSTNI 18)  *  EP8ARP 1 1  8 ) > I SGBARNI 1 31 /E )  0891 

TEFSTNI 18)  *  TEFSTNI 18)  ♦  EPCNKI 18)  0892 

OELEPNI 18)  *  0.0  0893 

GO  TO  1424,4261,324  *  0894 

426  OELEPNI 18)  *  EPCNKI 18)  -  EPCNPI 18)  0895 

424  CONTINUE  0896 

GO  TO  181  0397 

181  CONTINUE  0898 

IF! KLU2+KLU8) 266,266,267  _  0899 

266  CONTINUE  *  0900 
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LABEL  «  EFF. PLASTIC  STRAIN  0901 

FA  *  252626334743  0902 

FB  *  216263312360  0903 

FC  *  626351213145  0904 

WRITE  TAPE  8»KLU9,  NC,FA,Fd,FC,(EPBARN( J l 1 ,Jl=l,NC)  0905 

267  CONTINUE  0906 

IF(KLU2)25d,258,259  0907 

256  CUNTINUE  0908 

LABEL  *  TOTAL  EFF.  STRAIN  0909 

FA  »  634663 21 436 Q  0910 

FB  »  252626336062  0911 

FC  «  635121314560  0912 

WRITE  TAPE  8.KLU9,  NC,i  a.Fd.Fu . IT EFSTNI J  1) , Jl= l ,NC)  0913 

259  CONTINUE  0914 

00  862  125*1, NC  0915 

062  TtFSTNI 125 1  a  0.0  0916 

CALCULATE  NOME  STRAIN  Change  -  MaIkIX  DELEPK  -  FRAME  S12E  lfa5  X  l  0917 

00  191  III  a  1,NC  0918 

TEMPA  a  OELEPNtl  111/  SGdARN(ill)  0919 

M3  a  3*m  0920 

M32  *  M3-2  0921 

M31  a  M3— 1  0922 

0ELFPKI  M32)aTEMPA*UU  ?3l  (  1 1 l 1  *6lbUK I M32) S*ALl2  12 11 1 1 J *S  1 GUM  M3  l  0923 

1))  0924 

0eLEPK(M31)*TEMPA*(ALi22J(Ili)*iiouiUM31»-.5*AL1212(lll>*SIGUK(MV  0925 
HI  0926 

OELEPKI M3 J  *TEMPA*(3.*ALEA44Uill*alGUK(M3l )  0927 

191  CUNTINUE  0928 

CALCULATE  NODE  PLASTIC  STkaIN  -  MATRIX  EPSUK  -  FRAME  SUE  165  X  1  0929 

CALCULATE  NUDE  PUINT  STRAINS  0930 

00  192  123*1, NA  0931 

192  EPSUK  (1231  *  EPSUK  (I23J  ♦  UCLEPiU  1231  0932 

IF(KLU2+KLUd)260,260,26l  0933 

260  CONTINUE  0934 

LABEL  *  EFF. STRAIN  CHANGE^  0935 

FA  a  252626336263  0936 

FB  *  512131456023  0937 

FC  a  302145272562  0938 

WRITE  TAPE  8, KLU9,  NC,Fm,FB,FG, 10ELEPN( Jl I , Jl*l ,NC)  0939 

LABEL  *  NUDE  STRAIN  CHANGE  09*0 

FA  a  454624256062  0941 

FB  a  635121314560  0942 

FC  *  233021452725  0943 

WRITE  TAPE  8.KLU9,  NA,F A, Fb ,FC, I DELEPK (Jl) , J 1=1  ,NA »  0944 

LABEL  a  NUDE  INELAS. STRAIN  *  0945 

FA  a  454624256031  0946 

FB  a  452543216233  0947 

FC  »  626351213145  09*8 

WRITE  TAPE  8.KLU9,  NA,Fa,Fd,FC, ItPSUK  ( J 1 1  ,  J1  =  1,NA I  0949 

261  CONTINUE  0950 

C  CALCULATE  TUTAL  NODE  STRAINS  -  MATRIX  TOTEPS  -  FRAME  SIZE  165  X  1  0951 

00  201  U4al,NC  0952 

M32  «  36114-2  0953 

M31  a  3*114-1  0954 

M3  «  3*114  0955 

T0TEPSIM32I-EPSUK  tM32)+5IGUK(M52I/t  -GNU*SlGUK(M3l J/E  0956 


no 


w  a  a  o 


T0TEPSIM3D-EPSUK  t  M31  rlJl  )/t  -GNU*SItiUK(Mi2)/E  0957 

201  T0TEPS!M3  )*EPSUK  <M3  )*S16UMM3  1/iHRMOD  0958 

1FIKLU2 1 262.262, 263  0959 

262  CONTINUE  0960 

LABEL  =  TOT.  NUOE  STRAIN*  0961 

FA  =>  634663336045  0962 

F8  *  462425606263  0963 

FC  *  512131456260  0964 

WRITE  TAPE  8,KLU9,  NA,t-A,F8,f  u,  1  I0TEPSU1)  .  Jl=  1  ,NA »  0965 

263  CONTINUE  0966 

1FIKLU2 ) 268*268,269  0967 

268  CONTINUE  0968 

KLU9  -  3  0969 

KLUSI2  *  2  0970 

WRITE  TAPE  8, KLU9, KLUSIZ,  K,Pn,PM,PM,TlMfc  0971 

IFIN1NCLDI269, 269,273  0972 

273  CONTINUE  0973 

KLU9  *  4  0974 

KLUSI2  *  2  0975 

WRITE  TAPE  8.KLU9fKLUSI2iPrt.PM,Pri,^lNKCVtNINTOT  0976 

269  CONTINUE  0977 

IF(NINCLD>370,370,371  0978 

371  IFININCSW-2I422, 370,422  0979 

370  CONTINUE  0980 

IF  I  SENSE  SWITCH  6)421,422  0981 

422  CONTINUE  0982 

GO  TO  151  0983 

421  I F ( KLU2 1423,423,208  0984 

423  CONTINUE  0985 

205  KLU9  *  7  0986 

8  FA  777777777777  0987 

DO  206  125  =■  1,10  0988 

206  WRITE  TAPE  8,KLU9,KLU9,FA,FA*FA,FA,FA,FA, FA, FA, FA, FA  0989 

ENO  FILE  8  0990 

END  FILE  8  0991 

ENO  FILE  8  0992 

ENO  FILE  8  _  0993 

ENO  FILE  8  0994 

ENO  FILE  8  0995 

REWIND  8  0996 

IFISENSE  SWITCH  6)207,209  0997 

209  IF! J24-3) 210,207,210  0998 

207  CONTINUE  0999 

PAUSE  1  1000 

REWIND  11  1001 

CALL  CHAIN! 5,3)  1002 

210  CONTINUE  1003 

CALL  CHAIN  1 6,3)  1004 

208  KERRSW  =2  1005 

C  KERRSW  SET  TO  2  TO  MAKE  KLU2  *  0  AND  PRINT  A  CYCLE  1006 

GO  TO  151  1007 

998  KLU9  *  6  1008 

KLUSIZ  *  3  1009 

WRITE  TAPE  8, KLU9, KLUSIZ,  K, l8,PM,fcSUPRK,SGdARNt 181 , TIME  1010 

GO  TO  205  .  . . .  ......  .  .  1011 

ENO! 1,1, 0,0, 0,0, 1,1, 0,1, 0,6, 0,0,0) 


LISTS 


1013 


SUBROUTINE  OUTPUT! VALUE l, Si kPliNLUl) 

THIS  SUBROUTINE  SETS  KLUl  J  O  IF  IHc  CURRENT  CYCLt 
VALUE  =>  ABSFI  VALUED 
STEP  *  AB5FISTEP1) 

102  1F!VALUE-STEP)131,100,100 
100  NTESIl  *  IVALUE/STEPDl. 00001 
NTEST2  =  ( VALUE/STEPI+.SNa 
IFINTESTI-NTEST2) 131,130,111 

130  KLUl  *  0 
GO  TO  135 

131  KLUl  *  1 
135  RETURN 

ENOU, l, 0,0,0, 0,1, l, 0,1,0,0,0,0,01 


IS  TO  BE  PRINTED 


1014 

1015 

1016 
101/ 
101B 

1019 

1020 
1021 
1022 
1023 
1024, 
1025 


* 


CHAIN! 5,3) 
LISTS 


1027 

102B 
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C 451285  MATRIX  ANALYSIS  OF  INELASTIC  PLATE  -  LINK  5  -  WITH  CREEP _ 1030 

C  THIS  LINK  WRITES  A  SAVE  TAPE  EUR  RESTART  1031 

COMMON  TEFSTN,  TOTEPS  _ _ 1032 

COMMON  KLU4 •  KLU3,  KLU6,  KLU8,  NA,  NC ,  K  1033 

COMMON  KERRSW.  NINCLO,  KSET.  PM. _ Et___  GNU,  SHRMOO  _  _  1034 

COMMON  ALPHA,  BETA,  GAMMA,  TIME,  KLUISO  1035 

COMMON  PHTH,  KPMTM,  SIM.  S1J»  TS1GN,  TEPSN _ 1036 

COMMON  TALF12,  TALE23,  T ALE 31,  TALF44  1037 

COMMON  AL1212,  AL1223,  AL1231,  ALFA44 _ 1038 

0 1 MENS I ON  TALF12Illi,TALF2i(ll)( TALF311 11 ) .TALF441 11)  1039 

DIMENSION  AL1212<S5>,AL122J155),AL1231(55>,ALFA44I55I _ _1040_ 

COMMON  SIGUK ,  EPSUK  1041 

COMMON  SG8ARN .  SGBARP,  SGbARM.  EPflARN.  EPBARP,  OELEPN _ 1042 

COMMON  EPCNK,  EPCNP  1043 

EQUIVALENCE  I  TIMK1.TEFSTNI .  I  TOTEPS. OELEPK I _ _  1044 

C  THESE  FOUR  ARRAYS  ARE  EQUIVALENCES  TO  SAVE  CORE  SPACE.  1045 

!  DIMENSION  T 1MK1I 551. TEFSTNI 551. TOTEPS! 165 ) .OELEPKI 165 >  1046 

DIMENSION  PMTMI 10, 31, KPMTMUUi, SIMI165) ,SI4»  165*165), TSIGNI 11)  1047 

DIMENSION  TE.PSNI  1 1  l  iSIGUKt  1651  ,6PSUK(  165)  ,SGBAftN( 551, SGBARP 155) _ 1048 

DIMENSION  SGBARMI  55 ) ,  EPBARN1 55)  ,  EPBARP  1 55 1  yOELEPN  (551,  EPCNK  (551  1049 

DIMENSION  EPCNP155I  _ _ _ _  _ _ _  1050 

REWIND  11  1051 

WRITE  OUTPUT  TAPE  11.33 _ _  .  ...  _  _ 1052 

ENO  FILE  11  1053 

.  WRITE  TAPE  ll.KLU4,KLU5,KLU6.KLU8.NA.NC,KtKERRSW.NINCLD,KSEl,gft,Ei _ 1954 . 

1GNU, SHRMOO, ALPHA, BETA, GAMMA, TIME, KLUISO  1055 

WRITE  TAPE  iliPMTH  _  _ _  _  1056 

WRITE  TAPE  11, KPMTM  1057 

WRITE  TAPE  11, SIM  1058 

WRITE  TAPE  ll.SIJ  '  1059 

WRITE  TAPE  ll.TSIGN _  • _ 1060 

WRITE  TAPE  11, TEPSN  1061 

WRITE  TAPE  11.TALF12 _ _  _  _  _ _ 1062 

WRITE  TAPE  11.TALF23  1063 

WRITE  TAPE  11.TALF31  _  _  _  1064 

WRITE  TAPE  U.TALF44  1065 

_ MR1TE.TARE  U.*M212 _ 100& 

WRITE  TAPE  11.AL1223  1067 


MR  I TE  TAPE  11.AL1231 
WRITE  TAPE  11.ALFA44 
WRITE  TAPE  ll.SIGUK 
WRITE  TAPE  11, EPSUK 
WRITE  TAPE  ll.SGBARN 


1 068 

1069 

1070 


1071 

1072 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 


TAPE 

TAPE 

TAPE 

TAPE, 

TAPE 

TAPE 


11, SGBARP 
ll.SGBARN 


ll.EPBARN 
11. EPBARP 
11, OELEPN 
11. EPCNK 


1073 

1074 


1075 

1076 

1077 

1078 


WRITE  TAPE  11, EPCNP 
END  FILE  11 


ENO  FILE  11 
ENO  FILE  11 


SUBROUTINE  RUN  REWINDS  AND  UNLOADS  THE  DESIGNATED  TAPE 


1079 

1080 
1081 
1082 


1083 


PRINT  32 

1085 

KERRSW  -  5 

CALL  CHAINI6.3I 

1086 

1087 

32  F0RMATI//S2H  •  *  *  SAVE  TAPE  A-6  FUR  RESTART  AT  THIS 

33  F ORM AT  1 67HC0N T I N  99  SAVE  THIS  TAPE  FOR  RERUN  DECK 

POINT  *  *  *1 
45128  INELAST 

1088 

1089 

UC  ANALYSIS  ) 

ENOC liltOiOiOiOtif itOf I«0t0*0f0«0) 

1090 

*  CHAINI6.3I 

*  _ mu - 


1092 

^093 


"U3 


4S12tto  GRUMMAN  A I KC  RAF  I  ENGiNtERIrtu  CukP.  DtCK  NO.  45128  LINK  6 


C451286  MATRIX  ANALYSIS  OF  H«.4.*s)id  PlAIE  -  LINK  6  -  WITH  CREEP  1095 
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APPENDIX  VI 


STRESS  DISTRIBUTIONS  DUE  TO  UNIT  INITIAL  STRAINS 


The  basic  matrix  utilized  in  the  nonlinear  and  inelastic  analysis 
described  in  this  report  is  the  initial  strain  influence  coefficient 
matrix  [T  ].  Elements  of  this  matrix  provide  the  u^  stress  component 

“til 

caused  by  a  unit  initial  strain  at  the  v  stress  location.  The  first 
description  and  derivation  of  this  useful  matrix  vas  made  by  Denke  in 
195**,  Reference  20.  The  matrix  is  generated  by  a  structural  analysis 
which  may  be  of  the  force  or  stiffness  type. 

A.  Force  Method  Application 

The  method  utilized  here  is  a  simple  extension  of  redundant 
structure  analysis.  Essentially,  the  additional  vork  involved  is  the 
calculation  of  displacements  at  the  applied  loads  and  redundants  in  the 
statically  determinant  structure  caused  by  initial  strains.  These  dis¬ 
placements  are  combined  with  those  caused  by  the  applied  loads  and  re¬ 
dundants.  The  final  step  is  to  adjust  the  magnitude  of  the  redundants 
to  eliminate  the  total  displacements  at  the  redundants.  Fran  this  point 
the  determination  of  stress  distributions  for  the  redundant  structure 
is  carried  out  as  before. 

The  displacement  in  the  statically  determinate  structure  at  the  tth 
applied  load  or  redundant  due  to  initial  strain  can  be  expressed  by  use 
of  the  principle  of  virtual  work  as  follows: 


“tv  =  E/'xt-xr^ 

Volume  of 
ath  member 


(A-l) 


In  the  above  expression,  or,  is  the  displacement  at  the  tth  unit  applied 

load  or  redundant  due  to  the  v  initial  strain,  a  .  is  the  direct  stress 
^  Xtt 

in  the  a  member  due  to  the  unit  load,  and  c  is  the  initial  strain  in 
th  *v  v* 

the  a  member  at  the  corresponding  stress  point.  The  summation 

a 

indicates  that  the  virtual  work  in  all  members  affected  by  the  Induced 
strain,  c  .  muBt  be  considered. 


Equation  A-l  is  written  for  the  uniaxial  direct  stress  condition 
that  is  commonly  assumed  to  exist  in  bar  members  as  pictured  in  Figure 
(l).  The  effect  of  a  shear  panel  has  been  omitted  in  the  derivation 
for  simplicity.  Terms  necessary  for  Inclusion  of  shear  may  be  derived 
in  e  similar  manner. 


115 


Suppose  that  the  a  member  of  a  structure  is  a  bar,  with  cross 
sectional  area  A,  linearly  varying  axial  load  o^.A,  and  linearly  varying 

initial  strain 


Figure  1 


Then 


°rt  =  lflt(l  +  ^7} 
‘xv=  elv(l  -7}  +  e2v7 


and  after  evaluation 


where 


At-w* 

Volume  of 
member 


~  [Ylt  Y2t] 


L 

a21 


L 

a 


L 


a 


=  L 


11 

a22 

L 

12 

a21 

The  L matrix  can  be  similarly  determined  for  other  types  of 
“j 


structural  elements.  For  rectangular  shear  panels,  of  dimensions  b  and 
h,  it  can  be  shewn  that  the  La  factor  for  shear  strain  is  bh, 

the  panel  flat  plate  area.  For  nonorthogonal  structure  (swept  panels) 

other  .  factors  may  be  developed.  The  sum  of  such  matrices  for  the 
ij 

entire  structure  is  designated  [L^]. 
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The  matrix  expression  for  the  displacements  at  all  of  the  applied 
loads  and  redundants  due  to  unit  initial  strain  is 

[v^-NHV] 

The  e  's  are  the  Initial  strains  corresponding  to  each  of  the  member 
v 

loads  q^. 

Utilizing  Equation  A-2,  the  redundants  are  evaluated  by 


[v 


Of. 


rv 


t « 

i  rs 
i  J 


m 


=  0 


the  solution  of  which  is 


W  *  [r» :  r=v]  jy 


where 

j^sm  I  ^svj  [“raj  [arm  {  “rvj 
The  member  loads  in  the  redundant  structure  become 

W  •  [r-  i 

where 

flm  j  Flv]  "  [vj»  i  Vi,]  [l£,Tf^| 

Member  stresses  are  obtained  by  pre-multiplying  member  loads  by 
the  reciprocal  values  of  appropriate  bar  areas  and  shin  gages. 


Digital  computer  programs  which  are  available  for  conventional 

force  method  analyses  may  be  used  to  determine  the  T  ,  T  matrices. 

urn  uv 

This  is  readily  accomplished  by  redefining  several  input  matrices. 

(l)  replace  [y^L  the  usual  unit  load  distribution  in  the  static¬ 
ally  determinate  structure  due  to  applied  loads  by 


(2i)  x  (m  +  i) 

the  unit  diagonal  matrix  I  has  as  many  elements  as  there  are 
member  loads  in  the  structure. 

(2)  replace  the  usual  unit  load  distribution  in  the  static¬ 

ally  determinate  structure  due  to  redundants  by 


2i  x  r 


(3)  replace  [or^],  the  member  flexibility  matrix  by 


21  x  21 


The  values  are  the  geometrical  factors  of  Equation  A-2. 


Straightforward 


matrix  operation  will  now  yield  a  load  distribution  matrix  for  the  re¬ 
dundant  structure  which  may  be  identified  as  follows: 


EH1?] 


The  upper  portion  [r^  i  V  are  the  values  defined  by  Equation  A-3, 
and  stresses  may  be  obtained  by  using  Equation  A-4. 
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B.  Direct  Stiffness  Method 


Consider  a  truss  shown  in  Figure 
"locked"  (prevented  frctn  displacing). 


(2a)  where  all  the  nodes 
If  a  strain 


are 
in 


e  is  induced 
o 

particular  member  it  will  produce 

a  stress  in  only  that  member  equal 

to  a  =  -  Ee  since  all  the  nodes 
o  o 

are  locked.  This  stress  requires 
node  forces  as  shown  in  Figure  (2b). 

The  negatives  of  these  node  forces 
are  applied  to  the  "unlocked"  structure 


Consider  an  element  of  the  total  structure,  let  this  element  be 
supported  in  a  statically  determinate  fashion.  The  total  strain  e  is 
then  written  as  the  sum  of  two  strains  namely  -  the  strains  due  to 

stresses,  and  g  -  a  set  of  induced  strains.  Thus 
o 


U)  =  («c}  +  U0]  B-l 

Rewriting  this, 

{«CT}  •-  {«}  -  t«Q)  B-2 


U9 


The  stresses  may  be  expressed  in  terms  of  the  strains  by  using  Hooke's 
law, 

[a}  =  [b]{€o}  B-3 

Thus 


{o}  =  [b]{e}  -  [b]{e0] 


B-4 


(e)  (the  total  strains)  may  be  expressed  in  terns  of  the  displacements 
of  the  nodes  that  affect  the  particular  member  in  question.  That  is 


(«}  =  [a]{6]  B-5 

It  should  be  noted  that  the  expression  contains  the  basic  assumption 
governing  the  behavior  of  the  element.  For  example,  for  a  triangle  or 
bar  element  it  will  contain  the  assumption  that  the  strains  are  constant 
throughout  the  element.  Substituting  this  into  Equation  B-4 


{a}  =  [b][a]{6]  -  [b](eo]  B-6a 

=  [sd1M  -  [b]U0]  B-6b 

Fran  virtual  work  it  can  be  shown  that  the  nodal  forces  £f } , 
associated  with  (6),  may  be  expressed  in  terms  of  the  stresses  by: 

(f)  =  /[af{a)dV  B-7 

Vol 

Substituting  B-6a  into  B-7  yields: 

{f)  =  fl af[b][a]{6}dV  -  /caf[b]{« JdV  B-8 

Vol  Vol 


For  simple  elements  such  as  bars  or  triangles  where  the  strains  are 

assumed  to  be  constant  over  the  element  none  of  the  matrices  within  the 

integrals  will  be  functions  of  the  coordinates  x,  y  or  z  and  hence  they 

are  independent  of  the  integration.  The  first  portion  of  (B-8)  yields 

the  standard  stiffness  matrix  [k]  while  the  second  portion  yields  the 

nodal  forces  due  to  induced  strains  f *_} . 

oJ 

(fj  *  DOU)  -  V[af [b]{,0)  B-9 
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In  this  expression,  V  is  the  volume  of  the  element.  Noting  that 
[a]T[b]  =  [[b][a]jT  =  [sdf  yields 

[*}  =  CicHe)  -  vCsd]T{Co}  B-10 

The  total  stiffness  matrix  for  the  entire  structure  is  obtained  by 
superposing  the  stiffnesses  for  the  individual  elements  to  yield: 

(F]  =  [K]U}  ~  [Sd]T[V]{eoT}  B-ll 

Where  {F}  represents  all  of  the  external  nodal  forces,  {a}  represents 
all  of  the  nodal  displacements,  {e^}  represents  a  column  of  all  the 

induced  strain  components  and  [V]  is  a  diagonal  matrix  containing  the 
values  of  the  volumes  of  the  various  elements.  (Note  that  for  a  triangle, 
where  three  strain  components  may  be  induced,  the  volume  will  appear 
three  times).  [S^]  is  the  total  stress  matrix  array  which  is  obtained 

by  proper  arrangement  of  the  individual  stress  matrices  [s^]. 

Applying  boundary  conditions  (via  matrix  [BC])  to  equation  (B-ll) 
and  introducing  the  applied  external  loads  {P},  with  nodal  distribution 
[L],  yields: 

[L]{P}  =  [BC]T[K][BC]{A'}  -  [BC]TESdf[V][eoT}  B-12 

which  may  be  solved  to  give  the  allowed  nodal  displacements  under  the 
boundary  conditions: 

{A1}  =  [K^CLHPJ  +  C^l^tBcfCS^CVDte^}  B-13 

where 

[K^]  =  [Bcf[K][BC] 

For  the  total  structure,  B-6b  may  be  written  as: 


fop}  =  tSd]{A}  -  B-14 

where  [B]  is  a  diagonal  block  of  the  elastic  relations  [b]  for  each 
member.  Substituting  (B-13)  into  (B-lh)  and  using  the  relation  {A}  = 
[BC]{A'}  yields: 
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{aT]  =  [Sdl  CBC  3  [K^]"1  DL] 

+  -  [Bjjfe^}  .  B-15 

which,  to  use  the  previous  notation,  may -be  written  as: 

{a  }  =  '  [r  ]{p  ]  +  [r  ]{e  3  b-i6 

1  uJ  L  umJ''  nr  L  uvJI'  vJ 

The  first  matrix,  [r^j^of  this  expression,  is  the  conventional  distri¬ 
bution  obtained  for  unit  applied  loads  ad  indicated  in  Section  A  of  this 
appendix. 
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Mast  aerospace  structural  materials  exhibit  some  degree  of  anisotropic  strain 
hardening.  A  measure  of  anisotropy  is  introduced  into  structural  elements  by  such 
fabrication  processes  as  drawing,  cold  rolling  and  extrusion.  During  the  past 
few  years,  several  methods  have  appeared  in  the  literature  for  introducing 
inelastic  isotropic  material  behavior  effects  into  existing  matrix  analysis 
routines.  A  review  is  presented  of  one  of  these  methods.  It  is  essentially  a 
step-by-step  calculation  procedure,  and  corresponds  to  the  flow  theory  of 
plasticity.  The  method  has  been  extended  to  include  the  effects  of  anisotropic 
material  and  la  .formulated  as  a  standard  initial  strain  influence  coefficient 
problem.  Several  analyses  of  an  aluminum  alloy  (2024-T4)  •'hear  lag  structure 
which  has  been  tested  previously  for  the  Air  Force  are  carried  out  first  assuming 
isotropic  material  properties  and  then  anisotropic  properties.  The  resulting 
correlation  between  test  results  and  that  predicted  by  isotropic  theory  is 
reasonably  good. 

An  analysis  of  a  110CF  aluminum  shear  lag  structure,  carried  out  by  the 
incremental  method,  gave  reasonably  good  agreement.  However,  the  anisotropic 
creep  capability  was  not  checked  for  want  of  test  data.  The  approach  is  a 
reasonably  good  phenomenological  model  of  a  complex  physical  problem.  The  digital 
computer  program  submitted  is  suited  for  inclusion  of  other  msterial  nonlinearity. 
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